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Reihendarstellung axialsymmetrischer Geschwindigkeits- 
korrekturen und Strahldeformationen bei Windkanilen 
» fir den kompressiblen Unterschallbereich 


Von RoLanp FuUCHSHUBER, Nancy, France!) 


Ubersicht 


PRT NCR Te EY Ree Cem RYE OTE Rete NAL 


: Behandelt man fiir den kompressiblen Unterschallbereich mit Hilfe der 
_ Prandtlschen Linearisierung Modellumstrémungen in Windkanilen, so sind die 
an den Windkanalgrenzen auftretenden Randbedingungen zu beriicksichtigen. 
Es sei hier vorausgesetzt, dass die dem Modell zugehérige Str6mung in 
- einiger Entfernung von ihm in guter Approximation derjenigen einer gleich- 
" massigen Quell-Senk-Kombination2), eines Dipols oder einer einzelnen Quelle 
in einer Parallelstromung entspricht. Dann lassen sich die Wirkungen der 
endlichen Strahlgrenzen beim Windkanal nach der von PRANDTL-GLAUERT 
entwickelten Stérungsmethode interpretieren, und zwar als eine Zusatzstro- 
“mung, deren Geschwindigkeitspotential bei Uberlagerung mit dem zum unbe- 
_grenzten Raum gehdrigen Geschwindigkeitspotential die an den Strahlgrenzen 
auftretenden Bedingungen erfiillt. 

Ausserdem lisst sich so beim Windkanal mit offener MeBstrecke (Freistrahl) 
die durch das Modell verursachte Deformation der dusseren Stromlinie be- 
-rechnen. 

Eine mathematische Darstellung dieser Zusatzgeschwindigkeitspotentiale 
‘sowie der daraus resultierenden Zusatzgeschwindigkeiten und Strahldeforma- 
‘tionen mit Hilfe von einseitigen Fourier-Transformationen findet sich bei 
_ FucHSHUBER [1]?). 

In der Folge werden fiir eine Anzahl dieser Integraldarstellungen Reihen- 
‘entwicklungen angegeben, welche die numerische Auswertung in den meisten 
Fallen wesentlich vereinfachen. 

-- Zusatzlich sind einige dieser Reihenentwicklungen rechnerisch unter Ver- 
-wendung eines Magnettrommel-Rechners vom Typ IBM-650 im «Institut 
‘Universitaire de Calcul Automatique» der Universitat Nancy behandelt worden 
'(Programmierung vom Verfasser). 


1) Université de Nancy, Ecole Nationale Supérieure d’Electricité et de Mécanique. 
2) Das ist die Konfiguration einer Quelle und einer Senke gleicher Starke auf der z-Achse, 


jeweils im Abstand / vom Koordinatennullpunkt. 
8) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 288. 
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1. Berechnung der Geschwindigkeitspotentiale und der zugeordneten | 
Geschwindigkeiten | 


Fir das Gesamtgeschwindigkeitspotential der Modellumstromung ergibt. 
sich an der festen Windkanalbegrenzung (0 = 1) bei Ersetzung des Modells; 
durch eine gleichmassige Quell-Senk-Kombination (Rankine-Oval) im rota-: 
tionssymmetrischen Fall folgender Ausdruck: | 


osK _ UV 1 ae De ee : iii 
es R? 9 we AB J p L(e) oa p) sin(¢ p) dp ie ( ) 


Hierbei bedeuten: 

€ = z/R und 9 =7/R reduzierte Zylinderkoordinaten, 

R = den Windkanalradius, 
= V1— M? den Prandtl-Faktor, (M = Machzahl der Strémung ohne Modell), | 
= 1/R den reduzierten Abstand zwischen Quelle und Senke, | 


= das Volumen des Modells, 
= C/B; A=A/p. 


Die Bezeichnung der auftretenden Bessel-Funktionen erfolgt stets nach_ 
ERDELYI [2]. 
Der zugeordnete Ausdruck fiir die Gesamtgeschwindigkeit u23* = 0@92* dz’ 


ergibt sich dann zu: : 
co 


Ue el i wee s 
are gy / Tpy Sn b) cost 2) ap). 2) 


B 

A 

U,= die Anblasgeschwindigkeit, 
wg 

¢ 


Benutzt man die nach ERDELY!I [2], Gl. 7. 15 (64), abgeleitete Beziehung: 


t 
san $425) 35 Y1,n) (PF oe ad) ’ (3) 


wobei 7,,, die der Grésse nach n-te positive Wurzel der Gleichung. /,(¢) = 0 
mit ¢> 0 bedeutet, so findet man mit Hilfe von OBERHETTINGER [3] nach | 
einigen Umformungen folgende Ergebnisse: 


x co 1 eee ee al bet ss 
(+S acer rosinh(y, 9] fir 0 <b <A, 


poe OS pea os i 
+ a a 
| reel ¥1,n Jo(V1, n) 


A A 


e-Yine sinh(yy, , | tir A<2¢ oom 


4) Siehe [1], Gleichung III.1.14. 
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14> e~%1,n4 cosh(y,,, 6)| fir <4 >0, 
neat Jo (“1, ) k 


A A 


ea SA) ey 1+ ee LA) ty — 
’ : R 22 p 2 Py Jo(v1, n) i | ee Sees 


Sa ; Ce ein aS eae i al 
| 2. Jo(71, n) e~%1,n° sinh(y,,,4)| fiir ¢ > A. 


| 


In der Figur 15) ist nach (5) w98* fiir verschiedene Parameterwerte M und 
A graphisch dargestellt. Erwartungsgemass erkennt man, dass die Rand- 


q Figur 1 


R3 
u2S* als Funktion von ¢ fiir die Parameterwerte 
ORS eS) 


0 
A = 0,5 mit M = 0,866 und M = 0,968 
A=1 mit M =0,866 und M=0. 


Randgesamtgeschwindigkeit 


5) Die angegebenen M-Werte stellen auf drei Stellen hinter dem Komma genaue Naherungen 
der bei der Rechnung verwendeten Werte dar. 
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gesamtgeschwindigkeit mit wachsendem M rasch ansteigt, wahrend sie mit! 
wachsendem A, das heisst mit wachsendem Schlankheitsgrad des Rankine-- 
Ovals, absinkt. 

In der Figur 2 sind nach [1], Abschnitt IX, fiir verschiedene Parameter-: 
werte A bei festem M die Randzusatzgeschwindigkeiten u2>* fiir den Wind-. 
kanal mit geschlossener MeBstrecke und u32S* fiir den Windkanal mit offener’ 


; 
; 
| 


. 


- 025 = 
-190 -125 -{00 -075 -050 -025 0 G25 050 O75 {00 {25 [58 
€ ——> 
Figur 2 
3 
Randzusatzgeschwindigkeiten = 7 uesk beim Windkanal mit geschlossener MeBstrecke (Kur- 
0 


R : 
v u°SK beim Windkanal mit offener MeBstrecke (Kurven 2) als Funktionen 


ven 1) und Bye 


von ¢ bei konstantem M = 0,866 fiir die Schlankheitsgrade 4 = 0,5 undA = 1. 


MeBstrecke graphisch dargestellt. Erwartungsgemiss zeigt sich, dass die ab- 
soluten Werte fiir die Randzusatzgeschwindigkeiten beim Windkanal mit offe- 
ner MeBstrecke niedriger sind als die entsprechenden Werte beim Windkanal 
mit geschlossener MeBstrecke. 

In den Figuren 3 und 4 sind die Randzusatzgeschwindigkeiten u@S* bzw. 
Une beim Windkanal mit geschlossener bzw. offener MeBstrecke sowie die 
Randgesamtgeschwindigkeit u&3* beim Windkanal mit geschlossener MeB- 
strecke fiir die Parameterwerte M=0, 2=2 und M= 0,968, A=0O45 
graphisch dargestellt. (Die Randgesamtgeschwindigkeit ures beim Wind- 
kanal mit offener MeBstrecke verschwindet stets.) 

Man erkennt deutlich wieder den bereits festgestellten Einfluss von M und 
A sowie die gréssere Empfindlichkeit des Windkanals mit geschlossener MeB- 
strecke gegeniiber dem Freistrahl. : 


-Kurven 1: 


Kurven 2: 


: 
; 


“Kurven 3: 


~~ |? 


——e ee ie ies ee 
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Figuren 3 und 4 


R3 
Randgesamtgeschwindigkeiten TV beim Windkanal mit geschlossener MeB- 
0 


strecke als Funktionen von ¢. 


R3 
Randzusatzgeschwindigkeiten UW uesk beim Windkanal mit geschlossener MeB- 


0 
strecke als Funktionen von ¢. 


R3 
Randzusatzgeschwindigkeiten UV 40k beim Windkanal mit offener MeBstrecke 
0 
als Funktionen von ¢. 
Figur 3: M=0; A=2. 


Figur 4: M = 0,968; A= 0,5. 
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Die zu (4) und (5) analogen Ergebnisse fiir die Ersetzung des Modells durcl 
einen Dipol ergeben sich mit Hilfe der Beziehung (3) zu: 


oP, = AY aw | roy int)" 
0 


ur, = Oe ame | ripy cost P) a9) 
0 


Yi,n - ah = 
39 0 ate 3 Jo(7’1,n) ‘ asf (7 


2. Berechnung der Strahldeformationen 


a) Darstellung durch unendliche Reihen. 
Benutzt man die nach ERDELYI [2], Gl. 7.15 (62), abgeleitete Beziehung: 


il 1 t? 
Awan) Gt TEES APH WT ae (8)) 


wobei 7», die der Grésse nach n-te positive Wurzel der Gleichung J,(¢) = 0) 
bedeutet, so lassen sich die Deformationen der adusseren Stromlinie beim Frei-: 
strahl in folgender Weise darstellen: 


QSK V 2 " 
A : : 
i Re 218 é P 1y(P) sin(4 p) cos(¢ p) dp’) 
aie ¥ DEL at sy /[t potas oe 2 
E 2, Yo, n Ji(Vo, n) a sin fir ¢ << AS | 
Re 2nAB 


2 —— >= y an Ay ) e~770,n fiir = 
ee | 0,” 1\) YO, n 


0° 1 ae aie % * z : 


Das hier auftretende Integral ist im wesentlichen von BERTRAM [4] tabuliert: 
und von LENz[5] mit Hilfe des Liebmannschen Maschenverfahrens fiir: 
A= A/B = 1/3; 1/2; 1; 2; 5 bei konstantem 4 als Funktion von ¢/2 A sowie fiir: 
A = A/6 = 0 [Dipol-Problem, siehe Gleichung (10)]; 1/2; 1; 2; 5 bei konstantem Bi 


8) Siehe [1], Gleichung IV.1.2. 
7) Siehe [1], Gleichung VI.1.7. 
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als Funktion von ¢/2 8 berechnet worden. 


ea 1 - 1 s Cig 1 2 1 . A 
air zap J ee ae are 


(@>0). (10) 


pry 


4 ; : f love) ate 9) 
= a 2m | 2 +/ p a sin(¢ p) “| 
E 0 


= m, 1 i uF Oy ey erent ( > 0) (11) 


UR? 22 aa 2A0Nt Sion 
m,: Quell-Starke. 


Beachtet man weiterhin die nach SNEDDON [6] abgeleitete Beziehung: 


ae fir A> |C| 


co 


“lim > H(p) sin(A p) cos(é p) dp =: % 40) fir A= (C| (12) 


M-—1 
; 0 


O fir Aer it; (AceO)., 
so ergeben sich in Verbindung mit (9) und (11) folgende Grenzwerte: 


SETS deg ond 


R? 4nA 
GR a a, V iL ie 
pode hen Rt Bad fir A= |C|, ue 


Oar ate: (AS 0): 


[Der entsprechende Ausdruck fiir das zugehérige ebene Problem lautet zum 


_ Vergleich: 
eee far 4 >" |e 


, 


ee 3 
imag ai ftir A= C1, (14) 


Os tieaaelci; A> 0). 
F: Flacheninhalt der Modellkontur; H: Windkanalhohe.] 


8) Diese Beziehung findet sich bei G. Scumitz, Rechnungen zur Kanalkorrektur bei geschlos- 
sener und offener MeBstrecke, Forschungsbericht Nr. 820074 der Luftforschungsanstalt Braun- 
- schweig (1944), sowie im wesentlichen auch in [5], Gleichung (79, 14). 

8) Strahldeformation bei Ersetzung des Modells durch eine einzige Quélle. Siehe [1], Gleichung 
"VIII. 1.2, sowie [5], Gleichung (86, 13). 

10) Siehe auch [5], Gleichung (79,15). 

11) Siehe [1], Gleichung VI.2.5. 


er Shae Tee Te 
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itisenth ShGeine eNom 


~ TR 2m | 
: Q aa m it - 15)) 
Hime (6]=) _ 2 fir 6-0, | 
0 fir 6 <20c 


b) Geschlossene Darstellung durch analytische Naherungsfunktionen. 
Fiir die Strahldeformationen in (9), (10) und (11) hat man bereits eine} 

Reihe von geschlossenen analytischen Naherungsformeln aufgestellt. 
So ergibt sich zum Beispiel nach LENz [7]: 


| 


SK Yo, + 
Ax ~y Sap NS tanh( oe i) 
Ranke 2 


wobei hiernach im wesentlichen Ax fir 4 = 4/B =0 [Dipolproblem, siehe: 
Gleichung (17)]; 1; 2; 5, bei konstantem f als Funktion von ¢/2 graphisch wie-: 
dergegeben wird. : . 

Dividiert man in (16) auf beiden Seiten durch 4 und fiihrt man dann auf’ 
beiden Seiten den Grenziibergang lim durch, so folgt wegen lim Ax = Aa! 
nach elementarer Rechnung: 


Po ae ee gf fede 
Aas RY OnB sech ( 7 é). (17) 


Diese Beziehung wurde im wesentlichen zuerst von GRay [8] angegeben, aller-- 
dings mit 79; = 2,64 (v9, = 2,4048). Sie findet sich weiterhin bei BERTRAM [4], . 
p. 497 (dortige Fussnote 3). 
Schliesslich ergibt sich, etwa nach Lenz [7], S. 126, Gl. (16) und Gl. (17), 
die Naherung: 
Q m 1 Pats | 
Aa y — UR _— F + tanh gee :)| : (18) 


GLASER [5] (S. 317, Abb. 234, Kurve II) gibt im wesentlichen den graphischen 


Verlauf von A& nach (18) als Funktion von ¢/2 wieder. 
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Résumé 


Dans les souffleries 4 veine fermée ou libre, il faut tenir compte des corrections 
de vitesse et des déformations de la ligne extérieure de courant, dues aux conditions 


aux limites. 
Pour l’écoulement compressible subsonique, ces corrections et déformations sont 


représentées aux limites, dans le cas d’une symétrie de révolution, par des déve- 


loppements en série. 
Certains de ces résultats, aprés traitement numérique sur une calculatrice 


électronique, sont représentés graphiquement. 


q (Eingegangen: 26. Oktober 1960.) 


Ein supraleitender Gleichspannungsverstarker 
| fiir Helium-Temperaturen 


Von Suso GyGax, Ziirich*) 


3 Einleitung 
Die Messung sehr kleiner Gleichspannungen bei Heliumtemperatur stosst 
auf etliche Schwierigkeiten. Man kann wohl mit einem Galvanometerverstarker 
' [1, 2, 3, 4]2) Spannungen bis zu 10-® Volt messen, doch stéren die in den Zu- 
leitungen auftretenden Thermospannungen betrachtlich. Durch Umpolen der 
- MeBspannung mit einem Polwandler nach TEMPLETON [5] kann eine Differenz- 
- messung ausgefiihrt und die Thermospannung eliminiert werden. Dieses Ver- 
fahren wird jedoch um so schlechter, je kleiner die zu messende Spannung 
wird. PrpparD und PuLLAN [6] messen die Spannung direkt mit einem supra- 
_ leitenden Galvanometer, dessen Konstruktion jedoch enorme Schwierigkeiten 


oe 


7 


1) Institut fiir Kalorische Apparate und Kaltetechnik, Eidg. T echnische Hochschule, Zirich. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 297. 
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bereitet. Die wohl eleganteste Methode hat TEMPLETON [7] gefunden, bei der 
die MeBspannung bereits bei Heliumtemperatur so weit verstarkt wird, dass die. 
Thermospannungen der Zuleitungen keine Rolle mehr spielen. TEMPLETON und 
auch DE VROOMEN und VAN BAARLE [8], welche das Verfahren weiterausgebaut 
haben, benutzen dazu einen supraleitenden Modulator, bei dem ein supra- 
leitender Draht im Messkreis mittels eines periodischen Magnetfeldes die Gleich- 
spannung zuerst zerhackt, worauf sie dann wechselstrommassig hinauftrans- 
formiert und mit den iiblichen elektronischen Methoden gemessen wird. Doch 
sind besondere Massnahmen notwendig, um den Einfluss der Streufelder vom 
Modulator zu eliminieren. Anstelle dieses magnetischen Verfahrens kann man 
auch einen Warmechopper verwenden, wie er von uns entwickelt wurde [9]. 
Dass es jedoch auch méglich ist, bei Heliumtemperaturen gleichspannungs- 
massig zu verstarken, soll die vorliegende Arbeit zeigen. 


Arbeitsweise und Aufbau 


Die prinzipielle Arbeitsweise dieses Verstarkers, der die Ubergangskurve 
eines stromdurchflossenen Supraleiters im longitudinalen Magnetfeld benutzt, 
wurde schon friiher mitgeteilt [10]. Solche Ubergangskurven an Drahten von 
0.3mm @ einer 8°%-Pb—92%-In-Legierung bei 4,00° K fiir verschiedene Be- 
lastungsstréme zeigt Figur 1. Man bringe nun durch ein konstantes Magnet- 


J-IA J=100mA 


H[A/om| ——— 


Figur 1 
Widerstandssprungkurven von 8 At% Pb in In im longitudinalen Magnetfeld bei 
verschiedenen Belastungsstrémen. Drahtdurchmesser: 0,3 mm. 
Aim eieresers Lineare Approximation des Zwischenzustandes R = « (i — H,). 


feld geeigneter Grésse den Supraleiter in die Mitte des Zwischenzustandes. Jede 
zusatzliche Feldanderung wird dann eine Widerstandsainderung bewirken, die 
gemessen werden kann. Diese kleinen Zusatzfelder sollen nun durch die zu 
messende Spannung erzeugt werden. 

Figur 2 zeigt die von uns benutzte Schaltung. Die MeBspannung liegt an 
der Primarseite eines Transformators, dessen Sekundirkreis vollig supraleitend 


il Ni 


eek 
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ist. Ein solcher Transformator arbeitet als Gleichstromtransformator: im Se- 


= 
- 


. 


kundarkreis wird ein Strom induziert, der proportional zum Primarstrom ist 
und wegen dem verschwindenden Widerstand nicht abklingt, solange dieser 
Primarstrom konstant bleibt. Durch geeignete Dimensionierung des Transfor- 


-mators kann man somit den aus der MeSspannung gespiesenen, sehr kleinen 


it 


Gleichstrom zu betrdachtlichen Werten proportional verstarken und in den 
Spulen L, zur Erzeugung von Zusatzfeldern verwenden. Diese beiden ungleich- 


_sinnig gewickelten Spulen liegen je iiber einer Halfte eines im Zwischen- 
- zustand befindlichen supraleitenden Drahtes. Der Zwischenzustand wird er- 


 zeugt durch das Feld, das vom Belastungsstrom des Supraleiters in der Spule 
L, herrithrt. Ein Galvanometer liegt so tiber dem Supraleiter, dass es keinen 
_ Ausschlag zeigt, solange die zwei Halften des Drahtes denselben Widerstand 
haben. Erst wenn durch die Zusatzfelder der eine Widerstand kleiner und 
damit der andere gleichzeitig grésser gemacht wird, liegt am Galvanometer 
eine Spannung. 


eS eee NY 


ie 


i. Ts ee eT 


ow Vere se ee 


eer tite ee ee 
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Figur 2 
Verstarkerschaltung. Der eingefasste Teil befindet sich im Helium. 


Spannungsverstarkung 


Fiir die Spannungsverstarkung bekommt man (siehe Figur 2): 


oV OV of, Of; 
Sia hme i EMD 


mit V als Ausgangsspannung. 
Der erste Faktor wird einfach, wenn man das ganze Zwischenzustands- 


gebiet durch einen linearen Anstieg approximiert. Da wir in der Mitte des 
Zwvischenzustandes arbeiten, darf man dies, wie Figur 1 zeigt, sicher tun. Fir 


| 
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Pa re | 
H,<H <H, schreiben wir also: R = « (H — H,) und somit fiir die Spannung | 
V =a(H —4H,) J,. Das in den Spulen L,/2 von J, erzeugte Zusatzfeld ist | 
H = 1/2 n, Jz mit nz = N,/l (Windungen/Langeneinheit). Also wird jetzt: | 


or = te Joe | 
Der zweite Faktor berechnet sich dann aus der Transformatorgleichung 
Ly. iP = (L, + Ls) ? 
(wo Las = Ais OR die Gegeninduktivitat und k der Kopplungsfaktor) zu: 


ecg es 1 ft 
OS; I,+Ls V2 ray Es 
Ly Li. 


Besitzt die Spannungsquelle U einen verschwindend kleinen Innenwiderstand — 
gegentiber dem Eingangswiderstand R,, so wird 


Cle” 


“15 gms <2 


und somit die Spannungsverstarkung: 


yk ak ns Js 
15, 1s 
Ale eS 

Um nun eine gute Verstirkung zu erhalten, wird man eine moglichst 
grosse Steilheit « der Sprungkurve anstreben und mit einem grossen Bela- 
stungsstrom /; arbeiten. Figur 1 zeigt, dass die beiden Bedingungen miteinan- 
dergehen: je grésser der Belastungsstrom, desto steiler die Sprungkurve. 
Benutzt man statt einer reinen Substanz eine Legierung, so kann man einen 
hohen Restwiderstand erwarten, was wiederum die Steilheit vergrossert. Zu- 
dem kann das kritische Feld durch Zulegieren geandert werden: ein hdheres 
kritisches Feld bedingt in unserer Anordnung einen héheren Belastungsstrom. 
Figur 3 zeigt die kritischen Feldkurven der von uns benutzten In-Pb-Legie- 
rungsreihe. Als «kritisches Feld» H, wurde hier der Schnittpunkt des linearen 
Teils des Widerstandsanstiegs mit der Feldachse bezeichnet (siehe auch 
Figur 1). 

Fiir die Dimensionierung der Induktivitaten erhalt man folgende Forde- 
rungen: L, > L, oder N, > N 2, das heisst die Primarwicklung des Transfor- 
mators soll sehr viel grisser sein als die Sekundarwicklung, und die Kopplung 
soll méglichst gut sein (k ~ 1). Man wird somit als Kern ein Transformator- 
eisen mit grosser Permeabilitat verwenden miissen. Das erleichtert auch die 
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-Erfiillung der Bedingung L,, > Ls, welche besagt, dass die Zusatzspulen gegen- 
_tiber Ly, eine kleine Induktivitat haben sollen. 

_ Ferner soll noch R,, der Eingangswiderstand des Verstarkers klein sein. 
Dem sind bei einer Spannungsmessung Grenzen gesetzt: der Eingangswider- 
_ stand sollte ungefahr 10*mal grésser sein als der Widerstand, iiber dem die 


4 Spannung gemessen wird. Immerhin wird man bei Thermo- und Hallspannungs- 
-messungen R, sehr klein halten kénnen (~ 10-* 2). 


2 oe os 
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Figur 3 
Kritische Feldkurven der Legierungsreihe In-Pb. 


Nullpunktstabilitat 


Der erreichbare Verstarkungsgrad hangt wesentlich davon ab, wie genau 
der Betriebspunkt im Zwischenzustandsgebiet zeitlich festgehalten werden 
kann. Solange er sich auf dem linearen Teil der Ubergangskurve bewegt, spielt 
das bei unserer Anordnung keine Rolle. Erst wenn lings des Drahtes Asym- 
metrien im Widerstand auftreten, wird sich der Nullpunkt unserer Anordnung 
_ verschieben. Solche Nullpunktschwankungen kénnen sich ergeben: 1. bei Ver- 
_wendung von harten Supraleitern, die, wie unsere Untersuchungen an Ta er- 
geben haben, im Zwischenzustand sehr unstabil sind (siehe auch Justi [11] 
und Barrp [12]); 2. wenn sich iiber dem Draht ein Temperaturgradient aus- 
- bildet ; 3. bei Anwesenheit von dusseren Stérfeldern und 4. durch Schwankun- 


_ gen des Belastungsstromes. 
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1. kann weitgehend umgangen werden durch Verwendung eines weichen | 
Supraleiters, zum Beispiel einer passenden, gut getemperten Legierung; 
2. scheint bei uns nicht aufzutreten, und 3. kann durch eine Bleiabschirmung 
vermieden werden. 4. Die Schwankungen des Belastungsstromes kénnen be- 
trachtlich herabgesetzt werden, wenn man, wie aus Figur 2 ersichtlich, parallel 
zu Draht und Feldspule einen passenden Widerstand schaltet. Man sieht leicht — 
ein, dass dieser kleiner oder mindestens von derselben Gréssenordnung sein soll 
wie R + R;. Wegen M 
4 
Re Ree, 


wird mit unserer linearen Approximation fiir das Zwischenzustandsgebiet und : 
speziell fiir den Betriebspunkt R = 1/2 R,: 


1 
asd a mee ECT: ee, 
Ry, Ry, 


Fiir kleine Schwankungen muss also in der Tat gelten: 


R,< Kop Et. 


Ansprechzeit 


Wegen des kleinen Eingangswiderstands und der relativ grossen Indukti- 
vitaten ist eine grosse Relaxationszeit zu erwarten. Man findet: 


is ele ie pee SY 
R, LL, Ri eee 


t= 


Fiir ideale Kopplung (k = 1) und L, > Lz, das heisst hohe Permeabilitat, wird: 


Um einen-grossen Verstirkungsgrad zu erhalten, soll das Verhiltnis L,/L, gross 
und der Eingangswiderstand R, klein sein. Diese Forderungen fiihren aber 
gleichzeitig zu einer langen Ansprechzeit. Ein grosser Verstarkungsgrad kann 
somit nur auf Kosten einer langen Ansprechzeit erreicht werden. Ein Vergleich 
dieser beiden Gréssen ergibt: 


ae a Vr eee oe 


v ak ns Jy ENE 


und im Fall idealer Kopplung: 
PoP ph Ls 
u ans Jy VLo/Ly 
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Praktische Ausfiihrung 


NY 


_ Da wir bei 4,00° K arbeiten wollen, benutzen wir als magnetempfindliches 


= 


Element eine Legierung aus 8 At% Pb in In. Fiir diese Legierungsreihe ist der 
Verlauf der kritischen Temperatur mit der Pb-Konzentration in Figur 4 auf- 


diet 


‘y 
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‘Go 


36 
; 
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rel] 2 4 6 8 10 
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Figur 4 


Verlauf der kritischen Temperatur in der Legierungsreihe In-Pb. 


Figur 5 
Verstirker auseinandergenommen. Von links nach rechts: Trafo, supraleitender Draht mit den 
3 Potentialsonden, Zusatzspulen L,, Feldspule L,. 


SR ee ee ee ee Ce ee ey Ee ee ee 


“bessere mechanische Eigenschaften und gegeniiber Ta den Vorteil, dass er 
sich loten lisst. Mittels Pressen durch eine Diise wurde ein Draht von 0,3 mm @ 
erhalten. 3 Potentialsonden aus diinnem Kupferdraht wurden mit je 1,5.cm 


_Abstand angelétet (siehe Figur 5). 
Direkt iiber dem Draht befinden sich die einlagigen Zusatzspulen L3, die 


: 


| aus Létzinndraht (50% Pb — 50% Sn) von 2mm @ bestehen und je 2cm 


: gezeichnet (siehe auch [13)). Dieser weiche Supraleiter hat gegentiber In 


*?™ 


’ 
‘| 
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lang sind. Auch diese Legierung ist supraleitend und besitzt bei 4,00° K ein 
kritisches Feld von zirka 500 Amp/cm. Uber diesen Spulen befindet sich eine 
Nb-Feldspule mit Innendurchmesser 7 mm, Drahtdurchmesser 0,3 mm, Spu- 
lenlange 5 cm und N// = 135 Wdg/cm. Anstelle von Nb kann man aber auch 
einen diinnen Létzinndraht nehmen, dessen kritisches Feld ja gentigend gross. 
ist, um die Spule supraleitend zu halten; zudem kann man dann leicht léten. 
Der Transformatorkern besteht aus Eisenblechen, das Windungsverhaltnis ist 
N,/N,. = 700/1. Der Sekundarkreis besteht véllig aus supraleitendem Létzinn- 
draht, die Primarwicklung aus Nb mit 0,3 mm @. In den Sekundarkreis haben 
wir einen einfachen Schalter eingebaut: Uber einem kurzen Stiick des supra- 
leitenden Drahtes wurde eine Heizwicklung angebracht. Dadurch kann man 
einen Widerstand in den Kreis schalten, und ein eventuell vorhandener Dauer- 
strom wird aussterben. Zudem besteht damit die Méglichkeit nachzuweisen, 
dass eine Ausgangsspannung wirklich durch das Einschalten eines MeBstromes 
im Eingangskreis entsteht, und nicht durch irgendwelche Magnetfelder. 

Den Verstarker wird man mit Vorteil als Nullinstrument verwenden. Man 
erreicht dies einfach, indem man eine Kompensationsspannung tiber einen. 
geeichten Widerstand in den Eingangskreis schaltet. 

Die von uns verwendeten Widerstande und die aus der Geometrie berech- 
neten Induktivitaten haben folgende Werte: 


R,=810° 0, R.=1,21099, R,=1,110*0, R14 10-2 
N,=700Wdg, N,=1 Wag, N, = 15 Wag, u~50, 
I,~310Hy, L,~510-*Hy, L,~710-Hy. 


Mit diesen Daten berechnen sich die einzelnen Faktoren des Verstarkungs- 
grades: 


xr = 2,8 10-2 (gemessen: 1,2 10-3 Q) 
2 
2S, _ 95402 ( lee 
oy, = > gemessen: 1, ) 
A 
Oy : 
ie = 1210804, 


Damit bekommt man fiir 0V/0U = 8,4 108 anstelle des gemessenen Wertes 
von 1,6 10%. Die wirklichen Induktivitatswerte sind somit etwas verschieden. 
Die Nullpunktschwankungen am Ausgang wurden zu ~ 0,5 « V gemessen; 
dies ergibt auf den Eingang reduziert: ~ 10-® V. Die Relaxationszeit wars min, 
was der Berechnung entspricht. 
Zusammenfassend kann gesagt werden, dass man einen einfachen, kleinen 
Gleichspannungsverstarker fiir Heliumtemperaturen bauen kann, der zusam- 
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‘men mit einem iiblichen Galvanometer Spannungen bis zu 10-®V misst, bei 
einer Relaxationszeit von $ min. Eine Steigerung der Verstarkung durch pas- 
sende Wahl der Induktivitaétswerte wiirde eine Vergrésserung der Relaxations- 
zeit bewirken und ist deshalb wohl kaum wiinschenswert. Nur eine Vergrésse- 
Tung der Steilheit der Ubergangskurve hat keine Auswirkung auf die Ansprech- 
zeit, doch miissen dann gleichzeitig die Nullpunktschwankungen herabgesetzt 
‘werden. Ferner ist zu beachten, dass dieser Typ nur bei einer bestimmten, 
festen Temperatur arbeitet und zudem noch gegen dussere Felder abgeschirmt 
‘werden muss. Will man die Temperatur andern kénnen, muss man den Ver- 
starker in ein separates, kleines Dewar einbauen. 

_ Ich méchte Herrn Prof. Dr. P. GRASSMANN, Leiter des Instituts fiir Kalo- 
‘rische Apparate und Kaltetechnik an der ETH, fiir sein standiges forderndes 
Interesse an dieser Arbeit danken. Mein Dank gilt auch Herrn PD Dr. J. L. 
“OtsEN fiir die fruchtbaren Diskussionen und Ratschlage. 
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Summary 


-  <Asimple DC Voltage Amplifier using the resistive transition curve of a super- 
‘conducting wire of 8% Pb in In is described. This wire is held in the intermediate 
‘state by a constant magnetic field. The small current from the measuring voltage 
is amplified by a superconducting transformer and used to build up a small addi- 
‘tional field over the wire which in turn changes its resistance. This change in 
‘resistance is measured by a conventional potential drop technique with an ordinary 
‘galvanometer. The sensitivity of that device is 10-® Volt. It is limited by zero point 
fluctuations caused by instabilities of the intermediate state. 


‘(Eingegangen: 14, November 1960.) 
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Propriétés de convexité du type de Wey 
pour des problémes de vibration ou d’équilibre 


Par Jos—EPH HERsScH, Genéve?) 
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§ 1. Equations de Sturm-Liouville 


1.1. Nous considérons, dans un intervalle a < x < 6, l’équation 
L[{u} = (p(x)u’)’ — g(x) w+ Ao(x) u=0, (1) 


ot les fonctions g(x) et e(x) sont continues par morceaux, p(x) est Seep et 


est dérivable par morceaux, p(x) > 0, e(x) > 0; la masse totale Sot o(x) dx est 
supposée positive; et les conditions aux limites 


u'(a) —k, ula) =0 et u'(b) +k, u(d) =0; O<k,, ky <co. (2} 
Appelons 4, < A, < As < +++ les valeurs propres (carrés des fréquences 
propres: A, = w;); nous chercherons a décrire leur comportement lorsque l’on 
modifie quelques-unes des fonctions ~,¢,@ ou les constantes k,, k,. — Les 


fonctions propres correspondantes seront désignées par 14(x), u(x), t#9(x), ... 
on sait que (x) ne change pas de signe dans l’intervalle (a, b), prenons done 
U(x) > 0. | 


1) Institut Battelle, Genéve, et Ecole Polytechnique Fédérale, Zurich. 
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1.2. Le principe de Rayleigh 
La premiere valeur propre 4, peut étre caractérisée par le principe extrémal: 


A, = Min, R[v], ot 


b 
J (pv + q v2) dx + pla) ky v%(a) + p(b) Ry v*(0) (3) 
Riv] = 


b 
[evrdx 
a 


est le quotient de RAYLEIGH; on admet a concurrence les fonctions v(x) con- 
tinues et lisses par morceaux dans le segment [a, 0]. 
Une démonstration simple: 


(Puy 9% ye py Om) 


Uy Uy 


pv®+qv—Aov= puvr+quve+ 
= Puy 2)’ 2 Uy? oo Uy af Pie" ve ithe ne 
io v) + pe Legare. ss gr il ‘eon Baas 


d’ou en intégrant: 


b 
fot qe pov) ded [P18 oF] = — pO) By 270) — Bla) BA); 


a 


Q [p, g, — Ar; Pa, Ry; 41 | 
: (3’) 
= [ (v9 +908 — A008) de + Pla) be Ha) + P10) by (8) 2 0. 


a 


fon a l’égalité si v = u,, d’ou (3). 
Valeurs propres supérieures: A, est caractérisée par le minimum de R{v] 


sous les m — 1 conditions d’orthogonalité 
b 
[ouvar= pM 1 = (cf. (51),.p. 346). 


1.3. Variation de la densité de masse o(x): convexité de Ae 


Maintenons fixes (x), q(x), k, et &,; considérons notre probléme (1) (2) 
avec une fonction @,(x): valeurs propres Aled < Aled < As &:) < +++; et avec une 
fonction 0,(x): valeurs propres Ae) < led < en) < <.... Posons maintenant 
le probléme avec 1+ @2; § soit Ale+ oa sa premiére alent arog et u(x) la fonction 


2) Les chiffres entre crochets renvoient a la Bibliographie, page 321. 
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propre correspondante. u(x) est concurrente pour les problémes (03) et (Qe), done: 
Q (p,q, -— AP or; Ra ky; wl 20; | 
QP, 9, — AG? on; Ra kos ul 0; 
Q [p, g, — A * (0, + Os); Rar Ros ul =O. 


Multiplionsla premiére relation par A‘) /[A(@) + A{%)], laseconde par A) /[A‘) + Al], 
la troisiéme par. (— 1) et additionnons: 


Alert es) ese ee | “f (0, 2 Ph, u2dx>0; 
jien) a Ales) : 
L 1 : yr 
7 (0: +0) SS qe is aie) we SOILS CH) ILELATE: 
1 1 si i eee (4). 


Ae +02+ +++ +0n) = Aen) Ales) Alen) j 


Soit T la durée d’une période: A, = (2 2/T,)?, d’ou 


Tle: + es) ad Tle) aN Tle) : (4) 
(4) peut aussi s’écrire 
: : A u Wy 
(emey <2 (sien + Ger) (4") 
A 2 1 1 


ce qui signifie: Dans une famille linéaire 0(x) = @9(x) + x a(x), 1/A{?) = 1/A, (x) 
est une fonction convexe (vers le bas) du parameétre x). 

Remarque: Dans le cas particulier o(x) = const - @9(x), 1/A,(%) est une fonc- 
tion linéaire de x. 

L’inégalité (4) n’est pas nouvelle: cf. WEyL [11], p. 445; SourHWELL [10], | 
p. 483; CoLiatz [4], p. 314 («formule de DUNKERLEY»). — (4”) s’ajoute a des 
résultats de la théorie de Sturm, cf. [3], notamment pp. 57/58. 


1.3.1. Exemple: 
“u"+Ao(s)u=0, u(0)=—w(1)=0; 9, =1; a= cots (5) Daa 


(masse de Drrac). 
On obtient A) =m? et A) = 2atg (a/2); d’ot par (4): 


1 1 1 
ler + es) < 7 +3 


2 « tg (x/2) ’ 


3) Voir Post-scriptum a), page 320. 
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Ales +02) — gy 2 


voyons donc ce que vaut l’inégalité 


il if 1 
a2 «©=6 2a tg (aja) S n° 


Le membre de gauche vaut: 1/12 lorsque « = 0; (4 —2)/x? pour « = 2/2; 
1/x? pour « = x. 
1.4. Extension aux valeurs propres supérieures 


Soient 1, Ug, ..., U»_4,41 les premiéres fonctions propres pour (0, + @.); 
ans l’espace linéaire L(u,...., %,,4,41) engendré par ces fonctions, il existe 
au moins une combinaison w satisfaisant @ la fois aux m ponghor d’ortho- 
gonalité pour A\%), ef aux n conditions d’orthogonalité pour A{%),, soit en tout 
a (m + n) conditions. Nous avons donc a présent: 


a 
: Q [p, q; a (er 01; fs, ee u| = 0 ’ 
5 Q [p, q, — A), 09; ky; u| 20 , 
‘d’autre part, comme 4 € L(m,...,Unsn e son quotient de RAYLEIGH 


Riu) < Aes22, dot 
Op. 9g. — A, (or +00); har Bes ) < 


Multiplions la ates inégalité par Al), /[Ale, + (es) ), la seconde par 


Merten) ss a gles 
Si i = 1, on retrouve (4). 

L’adjonction dune répartition de masses Q2 augmente 1/A; (pour tout 1) au 
plus de la quantité 1/A\). 

; En raisonnant sur des équations intégrales, WEYL ((11], p. 445) obtenart des 
1912 Vinégalité correspondant a (5), ce qui contient en fait notre résultat, il Va 
utilisée pour l’étude du comportement asymptotique des valeurs propres. 


mie) (ate, + ae )], la troisiéme par (— 1) et tee 
roe ue, 
Beas = ae f (0, + @2) wdx >0; 

m+1 a ee 
: 1 1 1 
¢ Be. : Mm, n= 0} 12,3) vo. (5) 
1 pee I Ae, 
3 
Cas particulier n = 0: je poset = m +1; 
: ; £e ra ge ie (5’) 
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1.4. 1. Exemple: 
w+ Agu=0, u(0)=«(l)=0, | 


A 
Lb. “pout: O° <a <ee5 


iz 
1 
o1=45 pour fie ep eke Os 
al 
. pour a eS 


Ales + es) =; Jier* es) =472; Ale + es) =9n7?; Aes) =, Ales) 3 


(4) et (5) donnent 
Ae) <2 Aero) — 20%; AY <2 Aer = 18 22; 


tandis qu’en réalité 
Ae) ~ 16,47; Ae? ~ 96,56 . 


Remarque: Cet exemple montre qu’il n’est pas vrai que 1/A +) soit inférieur 
ou égal a 1/A) + 1/A@), car ici 
fe se i) ype A 
4m 2s 56,36 "06,50 4808 


Autre exemple contraire: 0; = 6,,, Q2 = 6,3 alors AG) = Ae) = + 00. 


1.4.2. Généralisation immédiate: 


1 : ub i 1 


= see es 
aaa TRE == (0) (22) +. (en) 
Ae ad Any +1 Aina Annet 


1.4.3. Adjonction de masses ponctuelles. — Si 0, est formée de k masses 


ponctuelles quelconques, on a A\%), = Alt), = +++ = 00 (car { 02 ule) vdx = 0, 


j=1,...,, entraine R@)[v] = oo). Posons, dans (5), m = k et m+ 1 = 7, nous} 
obtenons 


He <A (<A), $= 12,3000. 


L’adjonction de k masses ponctuelles ne peut pas abaisser A, ,, au-dessous de’ 
Vancien A,. 

Cette propriété se démontre directement ainsi: Soient 1,, Us, ug, ... les. 
fonctions propres pour (0, + @,); dans l’espace linéaire L(w,...,u;,,), il 
existe au moins une combinaison wu satisfaisant aux 7 — 1 conditions d’ortho- 
gonalité pour A” e¢ s’annulant aux k points portant les masses de 0, soit en 
tout 7 + k — 1 conditions linéaires; | 


Hier) — Ru] a Reedy] = Alen es) P condi: 
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- Notre énoncé est voisin du suivant bien connu: L’adjonction de & con- 
raintes ponctuelles ne peut pas élever A; au-dessus de l’ancien A; , ,. 

- Remarque: Pour un systéme vibrant guelconque dans l’espace a ie dimen- 
sions, on remplacera k par 3k dans ces formules. En effet, on a alors A‘) p= OO 
Théoréme analogue: L’adjonction de k corps rigides (systéme a 6 k degrés de 
liberté) ne peut jamais abaisser A; , 4, au-dessous de l’ancien j,. 


1.5. Variation des fonctions p(x) et q(x) 


o 
a 

4 

_ Maintenons fixes o(x), k, et k,; considérons le probisme: initial (1) (2) avec 
. et g,(x): valeurs propres A@»%) < Aut) < Aleut) < +++; et avec p,(x) et 
M.(x): valeurs propres A») < Jot) < jipoeo < .... Posons maintenant le 
Bisiine avec py ae Do, 1 aE Ges soient Mer+ by % + 2) <. Ales + Pay + 42) <a [ips + Par d+ 42) 
Dan . ses valeurs propres et #1, Ug, U3, ... les fonctions propres correspondantes. — 


Dans l’espace linéaire L(u,, ..., %m4n41), il existe de nouveau une fonction 
u satisfaisant a la fois aux m conditions d’orthogonalité pour Abe) ef aux n 
conditions d’orthogonalité pour 4? 2). donc 

‘a 

3 2 [pi M — pete; ky ky; u\ 20 ’ 

é . 

3 Q [ba Jar — Pet 0; ha hy: ul 0; 

‘ Q [br + Pa 9 + Yo — Amntt Qs Ra By; wl < 

‘Additionnons les deux premieres inégalités et retranchons la troisiéme: 


b 
{Albi tbe ds + a2) te (eet 4. Alen aa) ui owdx>0; 
a 


1 a > Alena) 4 yore)! —- m,n = 0,1, 2,3, .... (6) 


Le cas particulier m = n = 0: 


Abit Pa +2) = Meu i) sie Ae» qo) (6’) 


ee Oe ee Ce eee 


: 
se trouve également dans le livre de SourHWELL [10], p. 478-480; cette inégalité 
peut aussi s’écrire 


Pith, Ut% 
. Pie) Ore 1 ) 
A 2 2 > = (A(Pr ts) Be Ales 4) : (6") 


Dans une famille linéaire p(x) = po(x) + %- C(x), 9(x) = Go(X) + %° n(x), la 
premicre valeur propre A”? = A,(x) est une fonction convexe vers le haut du 


parametre x). 


r 
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5 

Remarque: Si €(x)/po(x) = const = n(x)/¢o(x), Ay(~) est une fonction 

linéaire de x. 
Cas n = 0: posonst=m+1= 1, 2,3,... 


APs Pe t+ 4a) = Apu 1) Bit Ae» ta) : (6")) 
“| 


1.6. Variation des liaisons élastiques k, et k, aux extrémités 


Maintenons fixes p(x), q(x), o(x); alors 2,, = 4,,(k,, %,); un raisonnement du 
méme type fournit les inégalités 


Ri) Et Rl) RO) A Re?) 1 4 4 
eas 2 ? 2 z= Q (en #4 (RY, RM) ae Avo (ke), RW?))\; 


(7) | 
Wi ae OR ANS AS 


Le cas particulier m = n = 0 exprime la propriété: Dans une famille linéaire? 
ki, =Ragt %Cy Ry = Ryo + % Cy, la premiere valeur propre A,(k,, k,) est une? 
fonction convexe vers le haut du parametre x. 


1.7. Application a la symétrisation des fonctions données 


1.7.1. Symétrisation des masses: Supposons a = — 6 (lorigine O est au: 
milieu de l’intervalle) et £(— x) = p(x), g(— x) = g(x), k, = k,. — Une répartitioni 
de masses 0(x) étant donnée, définissons la répartition symétrisée 


(2) = = [e(x) + e(— 2). (8), 


Appliquons (4") et (5) aux répartitions g(x) et o(— x), nous obtenons 


A) > Ae); plus généralement: A@),,, > Ale). (9); 


Remarques: a) Cette inégalité peut étre démontrée directement ainsi: ul?) 
étant une fonction paire, 2? = R@[ul?)] = RO [u®)] > ae, 

b) La symétrisation des masses définie par (8) est différente de la «symétri-- 
sation alternée» suivante (POLYA; HERSCH [6]): 

Définissons les deux répartitions symétriques@,(x) = (| x|)eto,(x) = e(—|x]) 1 
ona . 


0=20[f,9¢,-A4@ 0; fk; ule] 
ze Q Lp, Ue ae Oa Ra ka; ule) is Q LP, T= Ao Og) Ras Bon ue] ; 


l'un des deux termes a droite est donc < 0, d’ou 


A) > min (Aa), Alea) ; (10) 
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cette évaluation de A) est en sens contraire de (9): notre symétrisation des 
Beoasses par moyenne arithmétique (8) éléve A,, tandis que la symétrisation 
_alternée l’abaisse. ; 
Note: On démontre de facon analogue l’inégalité plus générale 

Meenea = min (AND, 122) . (10) 
(10) et (10’) restent valables, si p(x) et g(x) ne sont pas paires, avec 


y, P(x)=P(—|* |), gale)=a(l |), 9,(e) =9(—|% |); 


et, sl hee + Ry» avec n= hi k, 3 Rag = Rug = ha 


> 
= 
tos 
| 
= 
8 


— Encore plus généralement, la symétrisation alternée compare entre elles 
les valeurs propres pour des intervalles de longueurs différentes. 

c) Ces deux symétrisations sont essentiellement différentes de celles de 
B. ScHwarz [9] et BEESACK [1]. 

d) Reprenons l’exemple 1. 4. 1: 


0=0; AY ~ 16,47; AL) ~ 96,56 ; 
p25 5 Ma 2e8: Pas; IP = 182° 


On vérifie que A) > A et AY’ > AY); mais attention: Ae) Sb Ase): 


1.7.2. Des considérations analogues sont possibles pour: 
a) une symétrisation des fonctions p(x) et g(x) lorsque o(x) est paire et k, = F,; 
B) une symétrisation des liaisons aux extrémités: ha = ky = (k, + R,)/2 
lorsque p(x), g(x) et e(x) sont des fonctions paires. 
Ces deux symétrisations ont pour effet d’augmenter la premiére valeur 


propre. 


1.8. Application: Evaluation par défaut de A, a l'aide de la fonction de Green 


1.8.1. La fonction de GREEN g(x, &) (a < x,  < b) satisfait aux conditions 


aux limites 


SRE ee RCE Ee EE EN Ce ee ey ee eT ee 


g.(a, &) — ka g(a) =9,  by(b, €) + &y B(b, $) = 0 
et a l’équation différentielle 


Poder a mop 


ee OP 8 ee. 


ot 6, est une mesure de Drrac au point §. 
g(x, €) est donc, pour € fixe, solution de l’équation différentielle L ,[g(x, €)] = 9 


_ (cf. (1)) avec g(x) = 6; et Ay = 1/g(E, &); elle satisfait aussi (2). 


ZAMP X1I/20 
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Pour chaque & entre a et b, nous avons l’inégalité 
il 
— — 7 9:3 Ry, ky; 
O< Q [2. q, e(&, é) Og; a "“b v(x)| 


— 


. ” NG OIG 
+ 


q 


= [ pv? + qt) dx — TO). + pla) a) + 0) hy 2%) 


g(&, &) 


avec l’égalité pour v(x) = c g(x, €). On peut ainsi évaluer g(—,&) par défaut 


(ci. [24 ps277). 
Considérons maintenant une répartition de masses o(x) > 0, sa premiére 


valeur propre A‘) et sa premiére fonction propre’ u(x); posons v = u(x) dans — 


l'inégalité ci-dessus, multiplions-la par o(&) g(&, €) et intégrons par rapport a &: 


O< 


a 


b , 
[buts gut) de+ pla) by wa) + 210) ky | [ele 6.8) a8 


b 
= [ ol@ eG, 8) dé-@ p97, - =——* ‘eke 
q &, &) dé 


Admettons que g(&, &) soit finie (> 0); nous soustrayons de 


ale go aut > Bes Rp: ie 


J 0(&) g(& & dé 


Végalité 
Q[p,9, Ao; kk; u] =0; 
nous obtenons 
b 
; . 
ea -| / e #@ e>0, 
/ (8) e(&, &) dg | 


a 


b 
a <| o(€) g(é, €) dé .4) (11) 


a 


On al’égalité si la masse 0 est concentrée en un point €, quel qu’il soit: 9 = c- Og; 
alors A?) = 1/cg(é, €); la fonction propre fondamentale est alors simplement 
u(x) = g(x, §). 

L'inégalité (11) n’est autre que la limite continue de (4). 


*) Voir Post-scriptum b), page 321. 


\ fol. XII, 1961 Propriétés de convexité du type de WEYL 307 


= 


1.8.2. Corde vibrante inhomogene a extrémités fixées: 


u" +Ao(x)u=0 avec u(a) = u(b) =0. 


i 
4 (x — a) b= §) 

a Goa) pour *<& s i 

— (o — x) (€—a) ; ¢6g=% uss zi 

B. (b — a) pour x >é& 

dot 

; S : : 

A qo Sb—a | 0-8) E-aolé) a. (12) 


a 


(J apprends que cette inégalité est déja connue de plusieurs mathématiciens, 
notamment SCHIFFER, SZEGO et Lax.) 


1.8.3. Barre vibrante inhomogene, encastrée aux extrémités: 


4 
4 wlY — Zo(x)u=0 avec u(a) = w'(a) = u(d) = u'(b) = 0. 
. 

(11) reste valable avec _ 8 (e- 4) 


6 = oe 


a 
ie} 
= 


b 


r 


1 ss 1 : ' 
0 ~ 3(b—aP Geert Sere (13) 


a 


§ 2. Membrane vibrante inhomogene 


2.1. Nous considérons le probléme aux valeurs propres d’une membrane 
couvrant un domaine G du plan, et élastiquement liée le long de sa frontiére I’. 


Au + Ao(x,y)u=0 dansG, oe + R(s) w= 0 sur I’, 


ou n désigne la normale extérieure. 
Au lieu de (3’), nous définissons Q par 


O[—Aolx, y); Als); v= i/ [erad? v — Ap(x, y) dA + f h(s) v2 ds , 
G r 


ee ee, ee ee ee ee a an ee Se —- .. 


ou dA est l’élément d’aire et ds 1’élément de longueur sur I’. Le principe de 
RAYLEIGH dit que Q [— 4,0; A(s); v] & 9, avec Végalité si v = ,(x, y). 

Les inégalités (4) et (5) restent valables, la démonstration se laissant transposer 
sans difficulté; de méme, 
¢ (e hls) + Ral) ) ‘i . ; 
,] 2 s> E pales bio) 4 gles MOD] , (7’) 


m+n+1 2 


—— ee? ee ee oe 
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| 


Si le domaine G est symétrique relativement 4 une droite, nous pouvons - 
également transposer le § 1.7; une «symétrisation des masses» analogue | 
va étre indiquée ci-dessous. — En revanche, les résultats du § 1.8 deviennent 
triviaux, car la fonction de GREEN g(&, €) = o9¢ ici. 
Exemple simple: Membrane carrée avec k(s) = const; les chiffres indiquent | 
masse spécifique (ici constante par morceaux): 


A cee ay 


20 


0 Z 


l'inégalité de gauche s’obtient par symétrisation alternée (§ 1.7.1, b), celle de 
droite a l’aide de (4") comme au § 1.7.1. 


2.2. Membrane circulaire inhomogéne: symétrisation des masses, 
inégalité de NEHARI. 


2.2.1. Utilisons les coordonnées polaires (r, 9); soient G le cercle r < R et 
I’ sa circonférence; on donne g(r, 6) et (6); A, sera la premiére valeur propre, 
u, la fonction propre fondamentale. 

Considérons la répartition «symétrisée » 


1 6 
- aa fol 6) dO 
6=0 


et la liaison élastique (moyenne) 
ane: 
a rE) dO 


on a maintenant A, et %,.—- Comme %, (7, 0) = m,(r) est indépendante de 6, nous 
avons en vertu du principe de RAYLEIGH 


a 
ill 


Q[-Ave; Als); H)=Ol- Ae; &: H]=0< O[—Ao; A); Gl 
d’oti 
deg has (14) | 


(Cette notation est équivalente a: 
ise Be”) [ay] we oT [w,] > Ay .) 


are. 2. CaS Ue: k = oo (membrane a contour fixé) et 0(r) est monotone 
décroissante (dans 0 <r < R). 


La fonction de BESSEL Jolio7/R) (jo =2,4048) est concurrente pour le 
principe de RAYLEIGH avec 6; on sait que cette fonction est également décrois- 
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7 (3) do )] fig bday 2 
<A < R?1T,(¥1)) <5 ar 


as 1 
OU 0 = > per if o dA (valeur moyenne). Nous avons donc: 
-.- G 


Si o(r) est monotone décroissante, la premiére valeur propre est au plus égale 
= : . s 4 
54 celle d’une membrane circulaire homogéne de méme masse totale. 


: 2.2.3. Tel est par exemple le cas si o(r, 8) est une fonction surharmonique; 
“nous retrouvons ici l’un des résultats de NEHARI [7]. 


§ 3. Equation de Schrédinger 


3.1. Dans un domaine G de l’espace, de frontiére /’, considérons |’équation 


de SCHRODINGER 


q 2m , 
4 Au + =P [E — V(x, y, z)]u=0 (15’) 
avec les conditions aux limites 

* 

3 ou 

q —+hP)u=0 (Per); 


(x,y,z) est le potentiel, E le niveau d’énergie. — Posons A = (2 m/h?) E et 
(x, y, 2) = (2 m/h?) V(x, y, z), nous avons 


Au + [A— W(x, y,z)])u=0. (15) 


oe ee 
SN 


La premiére valeur propre A, est alors caractérisée par le principe de 
RAYLEIGH 


1, = Min, R[v] = Min, -°- pe seth peas _ 


ee ee ee ee ee ee eee 


ou dr est l’élément de volume, dS I’élément de surface. — Nous pouvons donc 
ici définir 


£Q (W(x, y,2) —Ay; A(P); 4] 


= III {grad? v og [W (x, y, 2) _ 7 v?} dt + {pao dor 0: 
G T 


wm 


ai% 


See 


_Q = 0siv = u, (premiére fonction propre). 

Remarque: Si v.est complexe, il faut écrire | grad v |? au lieu de grad? v et 
=|} v |2 au lieu de v?. 

r 


f 
- 
r 


, 
, 


;.  . = 
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3.2. Variation du potentiel 


Maintenons fixe k(P) sur la frontiére J"; varions W(x, y, 2). 


Premier cas: 

W, (x,y, 2); As) ; yu») . AW) ; uv) : Aw) é us) : 
second cas: 

W(x, y, 2) al”, ul), Aa), ws)» Aira) , us ; 
troisiéme cas: 


Wey.2= 5 (Writ Wi A, ts der tes Ay, ites 
3.2.1. Niveau fondamental d’ énergie 
QO IW, - AY; ky wm) > 
Q[W.—Ay"; k; ml > 


OW A be eae On 


Soustrayons la troisiéme relation de la moyenne arithmétique des deux pre-- 


miéres, nous obtenons A, — (1/2) (Ai) + Ai") > 0, c’est-a-dire 


(=: cle Ts) 1 p : 
Ag at A aha (16), 
(2 +Vo ) i : , { 
es As = (EV) + Es) | (16’) 
Nous avons de nouveau une propriété de convexité: Sz le potentiel V(x, y, 2) 


est de la forme V(x, y, 2) = Vo(x, ¥, 2) + x V(x, y, 2), alors le niveau fondamental 
d’énergie E, = E,(x) est une fonction convexe vers le haut du parametre x. 
On montre de méme, plus généralement : 


Si Cy, Cg») Cy 2 O et + eo+s* +e, = 1, | 
alors EQViAAVe+ 4609) > oy BU) + oy BUY 4 toy BY. fd 


3.2.2. Niveaux supérieurs d’ énergie. Posons 


~ 


W(x, y, 2) = ¢ W,(x, y, 2) + (1 — c) W,(x, y, 2), O<e<1; 


soient M4, Ug, »-. , Wm 4+n41 les premiéres fonctions propres avec W; dans l’espace: 


linéaire L(u,, ..., 4,41) engendré par ces fonctions, il existe au moins une: 
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e 
combinaison « satisfaisant a4 la fois aux m conditions d’orthogonalité pour 
ait, et aux m conditions d’orthogonalité pour 4\”%), soit en tout a (m + n) 


é Begone. Nous avons donc: 
Q [W, —AM); k; ul >; 
Q[W,—AW): k; u]>o 


n+1? 


paths 


? 


oo 


part, comme #€ L(u,...,%mninzi1), Son quotient de RAYLEIGH 
We) < Aes dou 

O[W = Aiea et ky ul< 

+ Multiplions la premiere inégalité par c, la seconde par (1 — c), la troisiéme par 

( 1) et additionnons: 


a Te eS 
% 


(Ati hee Cae st (1 —=% C) aay fff dt> 0 , 


Nee ee AT yA (18) 
HEM 0-9 © BW), + (1—e) BY. (18) 
Plus généralement: 
Sh 04, Oxy tase Oy oe 0 Etre, Fc, ee: oye, a7 
alors Elal+ a ae TT Ae a og Et : 


3.2.3. Application: cas d'une sphére. a) Soit G une sphére de rayon R; symé- 
trisons a la fois le potentiel V(r, 0, p) et le coefficient k(P) sur la frontiére I”: 


V(r, 6,9) =V() = 7 i f V(r 6, g) sind d0 dg: 


0-0 g=0 
be aa (PRP) sin d0 dq. 
y fe 


Alors 


~ 


Poh. 


oe ee ae tae ee Te eC Sella 4 kills Smad as tin dl pee 


En effet, soit 7, la premiére fonction propre avec W(r) et k; a, = ,(7), donc 


<O[Wr, 9,9) —; 26,9); mI=O(Wi)—A; Bs mM); 


_ soustrayons |’égalité 0 = Q [W(r r) — Ai k; u,|: nous obtenons 


: 0<( i, — 4) |] ac, 


ce qui démontre notre inégalité. 


eh i ed ee 


ZAMP_ 


b) Supposons par exemple k = oo oubien W(r > R, 0, y) = ov, c’est-a-dire | 
u =0 poury> R. Si en outre W(r, 6, y) est, dans la sphére, une fonction _ 
harmomique, } 


sal Wir, 0, y) dt; | 
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Ill 


Wir) W(0), done W = const = ( 5 


d’ot 
~ fag A 
1) = W+IO—-W+ Ze, car @, = — sin (7): 
We sa ob sald If Wir, 0,9) dt > (k =cosurl’). (19) 
3 tk (ee 


(N.B.: Sila condition sur [" était 0u/On = 0 au lieu de u = 0, le terme z?/R? 
serait 4 remplacer par zéro; l’inégalité est alors une conséquence immédiate du | 
principe de RAYLEIGH, méme si le potentiel n’est pas harmonique: il suffit de | 
prendre une fonction d’essai constante.) . 

Plus généralement, si W(r, 0, p) est sous-harmonique, l’inégalité (19) reste 


valable. — En effet, W(r) est alors monotone croissante; u(r) = r-) sin (x 7/R) 
est monotone décroissante et nulle sur J", donc 


We i < IT [grad? u Oke W uy” i at 


Te dt 
< fet aty fie Brat we. 
gadis Riles 


§ 4. Plaque vibrante encastrée inhomogéne 


4.1. Dans un domaine G du plan, de contour J’, nous considérons |’équation 
différentielle 
AAu — dA o0(x, y) 4 = 0 


et les conditions aux limites 
ou 


“w= —=0 sur I’. 
on 


Les inégalités (4) et (5) restent valables, ainsi que, si G a un axe de symétrie, 
(9) avec (8). 


4.2. Evaluation a l'aide de la fonction de Green 


La fonction de GREEN g(P, S) doit, S étant fixe, satisfaire A l’équation | 
différentielle Ap Ap g(P, S) = ds (mesure de Drrac) et aux conditions aux limites — 
g = 0g/Onp = 0 pour Pel. 


Py YSIS 
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ay 


frye 


g(P, S) = (1/82) rps Nn rps + (fonction de P, réguliére au point S). Cf. [8], 
=p. 140: Notre e(P, S) =I(P, S)/8 2; g(S, S) = s§/162; Porya et SzEG6 ‘ap- 
pellent Ss le «rayon intérieur biharmonique» au point S. 

A g(P, S) est la fonction propre (fondamentale) avec la masse spécifique 6, 
(masse unité au point S); la valeur propre correspondante est 1/g(S, S), car 


5 1 

q Apap er. eS gy 5 ese S}=0.- 

-Donc, en vertu du principe de RAYLEIGH, 

: ns [[(4e2a4 (40)? aA 

: S ew ! 2(5 , 

- [[9s0(P) \ dAp wv ( ) 

soit 

1%(S) <g(S, S) |] (dv* a4 (20) 


pour toute fonction v admissible, satisfaisant aux conditions aux limites; on a 
Végalité si et seulement si v(P) = c g(P, S) (cf. [2], p. 277). 

_ En particulier, introduisons pour v la premiére fonction propre u(x, y) pour 
la masse spécifique (x, y) ; 


I (Au)? dA = 2) I nada. 
u2(S) < ¢(S, S) / (Au)? dA = 4® @(S, S) i owdA. (20') 


Multiplions par o(S) et intégrons dans G, nous obtenons 


Peles eel 
7 < fats a(S, S) dA = sex ff e(S) 3 44s: (11’) 


/ za 


Sie ak, al 


on a l’égalité chaque fois que la masse @ est concentrée en un point. 


4.3. Plaque circulaire (de rayon R) 


4.3.1. Définissons de nouveau la répartition de masses symétrisée 


n= fe 


comme au § 2.2.1, nous avons 
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4.3.2. Cf. [8], p. 141: ss = (R? — | O'S |9)/R; donc par ani ‘| 


1 i 1 | 
Re = Ay = ne aie Arar fale 9 oe | 


r=0 


la détermination de Ae pose un probléme aux valeurs propres a wne dimension 

a l’aide du principe de RAYLEIGH, on peut obtenir des bornes explicites a 

gauche de (21). | 
§ 5. Problémes du type de Poisson 


5.1. Equilibre d’une membrane chargée 


Les méthodes qui vont suivre s’appliquent a des espaces 4 un nombre 
quelconque de dimensions et a des opérateurs différentiels plus généraux que 
ceux que nous considérons ici. — Nous nous occupons d’abord du probléme de> 
POISSON 


Au = — e(x, y) dans un domaine G du plan, w = 0 sur sa frontiére I’. 
On s’intéresse a l’intégrale de DIRICHLET 


u) = i grad?u dA = I[ e udA. 
vs a 


G 


Le principe de Dirichlet dit que 
vdA)" 
D(u) = Max UY : 


v=0 sur D(v) 4 


5.2. Variation des pressions: convexité 


Premier cas: 01; solution = u), 
Second cas: 023 solution ut, = Up, 
Troisiéme cas: © = (0, + Q2)/2 ; solution “= up. 


u est concurrente pour les principes de DirIcHtetT relatifs A 0, et a Qy, d’ou 


'D li ee ety a) abel] Ohtady 


VD(m) V D(x) 


> Dea) + YD) | > ls = Diu ; 


| Pg an) <> WM) + VD) || 


Sees 
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ne 


Si Von consideére une famille linéaire de répartitions des pressions 


AAA 


0(%, ¥) = Qo(*, y) + x a(x, y), 


VD) = V If om, dA 


st une fonction convexe (vers le bas) du parameétre x. 


alors 


. 


Pea 


Se SN 


= Remarque: Dans le cas particulier o(x, y) = const: @9(%, y), VD(u,) est une 
fonction linéaire de x. 

: 

e 5.3. Equilibre d’une membrane circulaire chargée: symétrisation des pressions 
y Nous imposons de nouveau uw = 0 sur la circonférence J". Symétrisons les 
pressions: 

z i 22 

E or) =z, | ol, 8) 46; 

- 6=0 

soient u(r, 0) la solution avec g, et u(r) la solution avec 0; 

j ——— eudA 6udA = 

3 /Dm > 12% gee =D), 

5 yD) yD) 

; D(i) < D(u) . (23) 
2 "7° , , 

4 5.4. Equilibre d’une plaque encastrée chargée 

; AAu = 0(x, y) dans G, u = 0u/dn = 0 sur le contour J’. Principe de Dr- 
-RICHLET: ii ’ 

4 A ; | (|[evaa 

y (Au)? dA = | owdA = Max, adnno or r\ 77 |” 

; ! I | ff (4p aA 

§ L’inégalité de convexité (22) se transpose de facon évidente. 

3 


5.5. Evaluation de u? et |] (ane dA al aide de la fonction deGreen 


Cf. §4.2: ApApe(P, S)=6,, g = 9g/dn,p=9 pour Pel. 
Lorsque la pression vaut 6; (mesure de Dirac: force unité au point S), alors 
a solution est g(P, S) et le principe de Dirichlet donne de nouveau (comme le 
‘principe de RAYLEIGH au § 4.2): 


g(S, S) > ee pourvu que v= 0v/dn=0 sur 1’; (20) 
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: 

d’ot en particulier 
u2(S) < g(S, S) if (Au)? dA = ¢(S, S) I oud , (24), 


avec l’égalité si et seulement si 9 = c 6s; alors u(P) = ¢ g(P, S) (cf. [2], p. 277% 
Multiplions par g(S) l’inégalité «(S) <//g(S, 5) V/ JJ oud et intégrons dans 
G, nous obtenons 


5 


| ff ouad < ff o(Syels, 8) dg, | 
I (Au)? dA = Ife udA <{ // o(S) ya(s, 8) a4.) ; (25); 


avec légalité si et seulement si @ = force ponctuelle = c 67, alors u(P) =c g(P,T).. 
(24) et (25) donnent 


w(B) < g(B, B){ |] (5) Ve(S. S) ash ; (26) 

on a l’égalité si et seulement si 9 = c dg, alors u(B) = c g(B, B). 
Cas particulier: @ = c dr, alors u(B) = c g(B, T); on retrouve (cf. [2], p. 277)) 
g°(B, T) <g(B, B) g(T, T) . (27) 


5. 6. Equilibre d’une plaque circulaire encastrée, de rayon R 


5.6.1. Si nous symétrisons les pressions o(r, 6) comme au § 5.3, nous; 
obtenons 


| (awraa < |] (au aa . (23’) 
5.6.2. Comme aux §§ 4.2 et 4.3, nous avons g(S, S) = s2/16 a et 
é R-|OSE. 
Ca R 


donc par (25) et (26): 


A : = ue = 
|f (4upad = ff onda < ie | | (RA) rar | or, 6) a ; (28)| 
r=0 
R 


6=0 

22% ; 

Ri- (ose | + | 

ae(B) | a, | (R — 7?) rar | ol, 6) d0. (29) 
r=0 6=0 


Par exemple, sig = 09, u(r, 0) = u(r) = g(r, 0) = (1/82) {77 In (7/R) + (R?— r)/2}. 
(cf. § 4.2 et [8], p. 141), tandis que (29) donne la borne | u(r) | < (R? — 7") /16 7oRt 
I oudA = u(0) = g(0, 0) = R?/(16 x), ce qui est exactement la borne donnée: 
par (28). 
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ay 4 Saath 


= § 6. Interprétations physiques 
6.1. Variations des masses sur une membrane vibrante 


6.1.1. Dans l’équation différentielle, A et o(x, y) apparaissent seulement sous 
a forme de leur produit A 9; cela nous permet d’affirmer: 
— Si,a module d élasticité E = 1inchangé, la masse spécifique devient = i,,0(x, ¥), 
alors la membrane aura une n-1eme fréquence propre w, = 1. - La membrane de 
module ad élasticité = E/A,, et de masse spécifique o(x, y) a une n-1eme fréquence 
propre w, = 1. 


6.1.2. Premiére valeur propre.— Nous considérons deux membranes vibrantes, 
-couvrant deux exemplaires d’un méme domaine G, et ne différant l'une de 
-l’autre que par leurs masses: 

3 Premiére membrane: module d’élasticité E, masse spécifique 0,(x, y), valeur 
propre fondamentale A(@”, 

Seconde membrane: module d’élasticité E, masse spécifique 02(%; y), valeur 

propre fondamentale A‘). 

_ Nous considérons alors les deux membranes auxiliaires suivantes: 

E Troisiéme membrane: module d’élasticité E/A%), masse spécifique @,(%, y), 
valeur propre fondamentale w} = 1. 

Quatriéme membrane: module d’élasticité E/A{), masse spécifique Q(x, y), 

valeur propre fondamentale wj{ = 1. 

La troisiéme et la quatriéme membranes ont chacune la fréquence fonda- 

“mentale w, = 1, donc également le systéme formé par leur réunion. 

Superposons la troisiéme et la quatriéme membranes et « collons»-les 
ensemble point par point: ces contraintes supplémentaires ne peuvent qu élever 
; les valeurs propres du systéme-réunion (cf. [5], p. 354): nous obtenons donc une 
 Cinguiéme membrane: module d’élasticité E(1/A~? + 1/A (¢)), masse spéci- 
-fique 0, + Qe, fréquence fondamentale > 1. 

Par conséquent, la . 
-  Sixiéme membrane: module d’élasticité E, masse spécifique 0; + @g, a une 
valeur propre fondamentale A+) > (1/A@) + 1/A?)-4, soit 


1 ce 1 " a 
les +e) ni Ales) a Alea) 


c’est l’inégalité (4). 
6.1.3. Valeurs propres supérieures 


Premiére membrane: module d’élasticité E, masse spécifique @,(*, ¥), 


_m-ieme valeur propre A‘2). 
Seconde membrane: ‘madule d’élasticité E, masse spécifique 0.(x, y), m-iéme 


Q 
valeur propre A‘), 
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Troisiéme membrane: module d’élasticité E/A'%), masse spécifique @,(x, y),| 
m-ieme valeur propre 1. 

Quatriéme membrane: module d’élasticité E/A%), masse spécifique g(x, y)») 
n-iéme valeur propre 1. 

Le systéme formé par la réunion de la troisiéme et de la quatriéme mem- 
branes, a (m — 1) + (nw — 1) fréquences propres @,, We, ... , Wm4n-2< 1, tandis 
que SON Oni 2-1 = Omin = 1. 

Superposons la troisiéme et la quatriéme membranes et «collons»-les en 
chaque point, ce qui ne peut qu’élever toutes les fréquences propres w; du 
systeme-réunion (cf. [5], p. 354): nous obtenons: une membrane (sur G) de 
module d’élasticité E(1/A'¢) + 1/A'%)), de masse spécifique = 0, +@., et de 
(m + n — 1)-iéme valeur propre > 1. Par conséquent: la membrane de mo- 
dule d’élasticité E et de masse spécifique 0, + @,, a une (m+ n— 1)-iéme 


valeur propre A+) > (1/A'%) + 1/A'%))-1, ce qui exprime précisément (5). . 


6.1.4. Remarque: La symétrisation alternée (§1.7.1b) s’'interpréte simplement — 
par le fait qu’une coupure de la corde ne peut qu’abaisser les valeurs propres. 


6.2. Varration du potentiel dans lV équation de Schridinger 


Considérée pour deux dimensions seulement, l’équation 


E Au + [Ao(x, v) — W(x, y)]u=0 


est exactement celle d’une membrane (module d’élasticité E, masse spécifique 
o(x, y)), dont chaque élément d’aire dA (au point (x, y)) est en outre rappelé | 
vers sa position d’équilibre w= 0 par un ressort infinitésimal de coefficient 
élastique W(x, y) dA. — Appelons A, la n-iéme valeur propre avec o(x, y) = E=1. 
On voit immédiatement: 

La membrane homogéne de module d’élasticité 1, ican spécifique =i, et 
«ressorts intérieurs» W(x, y), a une n-iéme valeur propre w? = 1. 

Nous considérons de nouveau deux membranes vibrantes (sur deux do- 
maines congruents) : 

Premiére membrane: module d’élasticité 1, masse spécifique 1, «ressorts 
intérieurs» W,(x, y), m-iéme valeur propre AW), . 

Seconde membrane: module d’ ween i masse spécifique 1, «ressorts inté-_ 
rieurs» W(x, y), n-iéme valeur propre Aw 

Nous considérons, a leur place, les fies anes auxiliaires suivantes: 

Troistéme membrane: module d’ élasticité 1, masse spécifique A), ressorts 
intérieurs W,, m-iéme valeur propre w2, = 1. 

Quatriéme membrane: module d’ élasticité 1, masse spécifique A‘), ressorts 
intérieurs Wp, n-iéme valeur propre w2 = 1, 

Systéme { troisiéme + quatriéme membranes}: a) 


m+n-1 — Onin 7 I, 


A NREL 
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__ Superposons la troisiéme et la quatriéme membranes et «collons»-les en 
chaque point: ces contraintes supplémentaires ne peuvent qu’élever w,,4,-1 
mcf. [5], p. 354): 

__ Cinquiéme membrane: module d'élasticité 2, masse spécifique = A) 4. A'"), 
pressorts intérieurs W, + Wy, On4n-1 =] 1. 

_ Sixiéme membrane: module d’élasticité 1, masse spécifique = (A) + Ai) /2, 
‘Tessorts intérieurs (W, + W,)/2, Onin-1 = 1. 

_ Septiéme membrane: module d’élasticité 1, masse spécifique 1, ressorts 
‘intérieurs (W, + W,)/2, valeur propre AlW.+W2l?] > (aM) 4. 40%9))/2, ce qui 
illustre V'inégalité (18) (ici par exemple pour c = 1 /2), ou (16). 


6.3. Variation des pressions dans un probléme du type de Poisson 


_ Considérons, par exemple, le probléme de l’équilibre d’une membrane 
_chargée: 
E Au = — o0(x, y) dans G, «=O sur la frontiére I’. 


(e(x, y) = pression; E = module d’élasticité.) Energie potentielle pour la posi- 
tion d’équilibre: E D(w)/2— [/ 9udA = — E D(w)/2. 
_ Premiére membrane: module d’élasticité 1, pression 0, (*, y), énergie poten- 
‘tielle = — D(u,,)/2 = — D,/2. 
- Seconde membrane: module d’élasticité 1, pression Q(x, y), énergie poten- 
“tielle = — D(u,,)/2 = — D,/2. 
- Nous considérons a présent les membranes auxiliaires suivantes (en équilibre 
également): 
Troisiéeme membrane: module d’élasticité 1/2, pression @,/2, énergie poten- 
tielle = — D,/4. 
 Quatriéme membrane: module d’élasticité 1/2, pression 0,/2, énergie poten- 
_tielle = — D,/4. 
|  Cinquiéme membrane: module d’élasticité k,/2, pression @,/2, énergie poten- 
tielle = — D,/4h,. 
- - Sixiéme membrane: module d’élasticité k,/2, pression g,/2, énergie poten- 
‘tielle = — D,/4 Rg. 
: Superposons la cinquiéme et la sixiéme membranes et «collons»-les en 
chaque point. Avant collage, l’énergie potentielle d’équilibre du systeme 
‘[cinquiéme + sixiéme membranes] était égale a — (1/4) (D,/Rky + Dy/ke); apres 
collage, le systéme a moins de liberté pour réaliser le minimum de l énergie poten- 
tielle, donc ce minimum (correspondant a l’équilibre) devient plus élevé: 
| Septiéme membrane: module d’élasticité (k, + k,)/2, pression (9; + @2)/2, 
énergie potentielle > — (1/4) (D,/Ay + D,|/ks). 

Huitieme membrane: module d’élasticité 1, pression (Q, + @2)/2, énergie 
potentielle = — D(u, +0/2)/2 2 — (Ry + Re) (Dy/ky + Do/ke)/8. 
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Nous avons donc |’évaluation 


\ D D k k 
gee city (ky + ep) x a 7) = Dit Dek Dig Pan age 
2 


= (Dicey Dy +(V/D.5 - ry. 


Cette évaluation est la meilleure possible si l’on choisit k,/k, = /D,/D,, on 
retrouve alors bien l’inégalité de convexité (22). 


Post-scriptum 


a) G. Porya et M. ScHIFFER: Convexity of functionals by transplantation, 
J. d@’Analyse Math. 3 (2° partie), 245-345 (Sie) notamment p. 286-290, 


ont établi des propriétés de convexité pour Bs A; et pour Dy A; qui généralisent — 


(4), (6’) et (16) dans un sens différent de (5), (6) et (18), c’est-a-dire autrement 
que les inégalités de H. WEYL. 

Voir aussi dans cet ordre d’idées mes deux notes: Caractérisation variation- 
nelle d'une somme de valeurs propres consécutives; généralisation d’inégalités de 
Pélya-Schiffer et de Weyl, C. R. Acad. Sci. Paris 252, 1714 (1961); et Inégalités 
pour des valeurs propres consécutives de systémes vibrants inkomogeénes, allant 


«en sens inverse» de celles de Pélya-Schiffer et de Weyl, C. R. Acad. Sci. Paris 252, | 
2496 (1961). — Voici quelques exemples de ces inégalités: 

sia S ia a i < fee e Tar , (30) 

—— ne ba ROE Ss rr i: rar (31) 


(l'inégalité de gauche est de POLlyA-SCHIFFER). 
L’exemple considéré ici aux §§ 1.4.1 et 1.7.1 (d) illustre bien ces inégalités: 

(30) donne 
1 = 2 ze 1 1 


<5 <—=+ 
ne ie) 72 4 72 


soit 1,6. 1? < Ay < 22? (en réalité A) ~ 16,47); (31) donne 


1 1 
Pen eal SA Brion perme oy i 
1 2 


ot ; 
m0 4 7? 9 7 16 x? 


y 


i 
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soit 


~*~ 


-c est-a-dire 


5 eet og, 7205 
Sat aed ° ge) ~ 288 x2 ’ 


a t 1 
0,0633 SS lew -b ye) < 0,0721 
f 1 oe i) 1 


+ <. 


en réalité — ~ 5 
( ated Jie) 16,47 96,56 


~ oe 


Dans l’exemple du § 1.3.1, A: *%) = AM) = 4; (30) donne donc 
1 1 1 1 1 
ee ee Beare as ro 
Zoe = 
soit 
3 1 1 = 1 
2 ah te te) 
4n a 2atg oO I 


we 


Spour0<a<x. 


Za 


b) L’inégalité (11) peut étre précisée par la relation connue 
fo 2) 
= 


A 


t 


io — | eee 448, 


" qui est une conséquence du théoréme de MERCER, cf. R. Courant et D. HIL- 
BERT [5], et B. ScHwaRz: Bounds for sums of reciprocals of eigenvalues, Bull. 
Res. Counc. Israel, 8F, 91-102 (1959), notamment p. 94. Une interprétation 
| en est indiquée dans ma note citée en dernier lieu ci-dessus. — On peut ainsi 
" préciser également (12), (13), (11’), (21); on comprend pourquoi on ¥ a l’égalité 


si et seulement si toute la masse est concentrée en un point: alors A, = Ag == 
.++ = + oo, la somme infinie se réduit a son premier terme. 
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Summary 


General inequalities of type (5) for eigenvalues were established as early as 1912 
by H. Weyt [11] through consideration of integral equations; he used them for 
studying the asymptotic behaviour of 4,, when ” + oo. These inequalities seem to 
have been forgotten, as some of them were independently re-discovered by several _ 
authors (see [4], p. 314 and [10], pp. 483-484). — Such inequalities are here derived 
and applied in various ways to special problems (vibrating strings, rods, membranes 
and plates). Explicit lower bounds for A, are obtained in a somewhat similar manner 
in terms of the GREEN’s function. For ScHRODINGER’s equation (15’), the very 
simple inequalities (17’) are found. Analogous inequalities are established for a 
class of equilibrium problems. Many of the inequalities thus obtained express 
convexity properties (cf. the Post-scriptwm above and the paper by Pérya and 
SCHIFFER quoted there). 


(Regu: le 3 décembre 1960.) 


Green’s Function for Laplace’s Equation in a Circular Ring 
with Radiation Type Boundary Conditions 
By James F. Heypa, Cincinnati, Ohio, USA}) 


1. Introduction 


Let (7,0) be the polar coordinates of a point of the circular ring 
0 <7 <7 <-y%, and consider the problem of solving Porsson’s equation, 


YT = 0?T 1 oT i, STE 2 


oe toa ta oe oF), (1) 


1) Aircraft Nuclear Propulsion Department, General Electric Company. 
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‘subject to the boundary conditions 


mT) = bien A) (20), 2) 


oy the Ne, = gels, 0) (ha = 0), (3) 


where the source function F and the functions ¢,, ¢, are given. 

Equation (1) describes the flow of heat in a long hollow cylinder, the flow 
‘taking place in planes perpendicular to its axis, with ; = H;/k, where H; is 
the coefficient of surface heat transfer at r = 7; and k is the thermal conductivity 
of the cylinder material. We shall assume the /; constant. 

When ¢, = ¢, =0, the boundary conditions (2), (3) describe cooling by 
forced convection and are sometimes referred to as radiation boundary con- 

ditions. 

The solution of the problem defined by equations (1), (2) and (3) is 


DAP AOS | Ge. 0; 79, 99) F(%, 99) 4Ag + [bro O):G(r, 0; %, On asa, 14) 
R G 


"where G(r, 0; 79, 9) is the appropriate GREEN’s function; dA, is an element 
of area of the plane annulus R: 7, << 7 < 17; C is the complete boundary of R; 
the function ¢ = ¢; when 7, = 7; and ds, is an element of arclength on C. 


2. Derivation of G(r, 9; ro, 9%) 


It is convenient to denote the point (7, 6) by the vector z=7 e’® and to 
“write G(r, 0; 7, 9) = G(z, %). The GREEN'S function G(z, %), where 2, 2) are 
distinct points of R, is defined to be the solution of 


V?G(z, 2) = —0(2 — 2) (5) 


with the boundary conditions 


It may be written 
G = G* 1s = , (7) 


where the singular component S is given by 
S(z, %) = —z~ Re log (r e** — roe) (8) 


and the non-singular component G* is a solution of V?G = 0 subject to condi- 
tion (6). 
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One shows readily that 
1 cr 1 (%\" 
S(é, 20) = z— F logy ek (“") cosn (9 — a) (¥9 <7) 


il 


Fa [loam + 3% (Z) comm (0-09), o>, 


I 


while the most general G* in terms of circular harmonics is 
G*(z, 2) = Ag+ B, logy + Dd" (A, m+ Bir") cosn 0 | 
1 
+ D(C, 4+ Dr-*)sinn 6, | 
1 


m nm 


wherein the ‘constants’ Ag, By, A,, B,, C D, (n=1,2,...) are now to bes} 


determined so that 


a 1h, S| Ga te ee 


Finding and simplifying these ‘constants’, although somewhat of a pro- 
digious task, is well worth the effort in view ei the reasonably simple results. 
One finds, upon introducing the substitutions 


=e", Baem, Bien OSnmSm), (12) 
am y vy a 
a 
: (i, ines ~) : i. - a) | 
Ags = " J (13) 
hy hg No - = B+ — 
Ys 
1 
h is (hs =o) er a 
By = 5+ ——_—_____=) (14) 
Poviies rap h, 
1 2 Ng 2 + Yo 
h 1 n 
en ) sinh —(i =>) I | 
hip ES Nt td ae el ee (15) 
: HU AS » 
: (= - + fy hy) sina ny +n ( +28) cosh n . 
’ Yo 1; 
hy IN oe 1 n 
an Vi COS Dy hy ( n — 7 sinh (72 — mp) + a (1, — zh coshn (2 — m9) 


Zo n* hy 
("- + hy hs) sinhn 7H, + 2 (“4 + ) cosh” 1 
173 Yo V+ 


(10) 


D 
A Re tans, (17) 7] 


| 
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Ba alate 


_ The final form of the GREEN’s function G(r, 6; 7), 89) may now be obtained 
by substituting the above evaluations of Ao, By, A,, B,, C,, D, into equation 
(10) and combining the result with equation (9) in accordance with relation (7). 
After much algebraic refinement, the final result turns out to be: 


A 


a = . 
For eS Pe ee oe hes, tor “0°Ss Hy = 1) So 5 
ar 


G(n, 9; Mo, 9) = Zs | ‘ * “) ‘ ee a) 


h h 
hy hy Ny + = 9 a 


2 


- n : 1 
“4 co (= coshn my + hy sinhn no) (—- coshn (yn, — ») + os sinhx (7. — n)) 
1 2 


+72 es i 


: h 

; ae +h, Zs) sinhn ny + ” ise + T)eoshn Ne 

Since G(z, 25) is symmetric in z, 9, the form of G when 7, S17 < 7% S718 
_ obtained from equation (18) by interchanging 7 and %p. 

A 

3 3. Green’s Function Vanishing on the Ring Boundaries 

4 By letting h, > oo, hy > oo in formula (18) we obtain 

[ 

: ; yg) eis = py sinh nm  sinhn (12-1) oon (9 — 
BG(7, 9; No, 0) =Z— - + Pd, aes “ cosn (0 — Oo) |, 


a reasonably simple form of the GREEN’s function which vanishes on the outer 
circle 7 = ro, i.e., 7 = 2. Equation (19) applies when 0 S m% <7 <=7,; the 
applicable form when 0 < 4 < mp S mz is obtained from (19) by interchanging 
n and 7, thereby yielding the GREEN’s function which vanishes on the inner 
Ecircle 7 = 0. 
; This case is treated by HirBert and Courant [1]*) who derive the equi- 
valent of (19) using complex variable methods. Although their analysis is 
- short and painless, their final result, expressed in terms of theta functions, does 
- not lend itself easily to analytical computations. It is of some interest to show 
_ that the two results, superficially very dissimilar, are indeed equivalent. 
HitBert-Courant take 7, = g¥?, 7, = g-1?, (0 <q < 1), so that 47% = i 
Free of the theta function notation, their result is 


cosn (8 — Oy)). 


(18) 


(19) 


G2, 4) =—z2 E log + log | 27 '°8*0°8?| + log| 2 . ps og) 0] ie 


2) Numbers in brackets refer to References, page 327. 
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where 0° 
IT aj od ) (1 Ope a) 
fa Ge Zz 
Q SF a eS qe”) 
[[a-@ 44 (1 4 {—) 
st Zo 2 


It is sufficient to compare (19) and (20) for 2 real, i.e. 2 = 7%, 09 = O or Tes. 
To be specific we take 0, = 0. Equation (20) can then be rewritten, when) 
Py 0, aS 


nD fgogre loge rn] 
CaN ey [ oe ee 


: : (22) 
if Y 
— 5 log|1 + (7°) — 2 7 cos 6] — log |Q]. 
It is then a simple matter to identify the term 
(pL SAE Lace Na op a 
7a elas ae ert 
with the term 7 of equation (19). It remains then to show that 
Noe — 1/2 q 
3: sinhn 7 sinhn (nj, — 7) ; 
72, ss msinhn 7, api ad | (23) 


1 aa 
= Go log [1 a (2 y —2 (< 2) cos | logy e | 
On the right side of equation (23) we have, for the first term, the known 


expansion co 
~> > log [1 + (<2 y —2 “! cos 0] = a, ( “! ‘ve beac (24) 


n 


From equation (21), 


19g|0| Tog [J 1-7 2 + og [1 = v2 
v=] y=1 | 


ia =|. 


— tog [J [1 —A-* ry2 | log TT | 
y= 


Consider the first term of equation (25). We may write 


log [J |1— fv £| = Slog 


ve=1 | vel 


Be et rey ie N= 
yee ee Tea (re 


4y 2 
t= 7; ie 
0 
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Bet G ee tm OF PS 1 ey 


z vy \” cosn 0 
) / | — vty = y Y s Y An 
Is 1 % y eS % (+) ——S=— 
Se v=1 0 pat a= 0! n 


oo 
ey Ket (<.)" cosn 0 
1-4 \al oe” 


n=1 
where we have used an expansion as in (24), interchanged orders of summation 


g that 7, = — 2 log 71, we may write finally 


log [] 
v=) 


es Ine PCr: (27 


d nm sinhn 15 


Proceeding similarly with the other terms of (25) one easily obtains the 
esult 


‘6 
co 72" cosh (n log— 
, 


7” : 
¥ a ee cosh (” logy 79) 
; log | Q| wel ac cos” 0 po ne ea cosn @, (29) 


simplified, the result is the required left member of (23). 

Other specializations of the general form (18) are also of interest. Thus 
‘one might keep /, finite and let h, > oo, as would be required in determining 
‘the steady temperature in a long hollow cylinder with cooling at the inner 
‘surface and temperature 0 at the outer surface. 

The form (19) is only slightly more involved when the ring boundaries are 
‘non-concentric. This result, involving bi-polar coordinates, was previously 


presented by the writer in reference [2]. 


Si it is a simple matter to show that when (29) and (24) are added and 
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Zusammenfassung 


Es wird die Greensche Funktion der Laplaceschen Gleichung fiir einen Kreisring | 
hergeleitet, wobei Randbedingungen vom Strahlungstyp angenommen sind. In} 
demjenigen Spezialfall, in dem die Randbedingungen das Verschwinden der 
Greenschen Funktion verlangen, wird gezeigt, dass die gefundene Darstellung der 
Greenschen Funktion iibereinstimmt mit der Formel, wie sie im Buch von HILBERT- 
CouRANT unter Verwendung von Thetafunktionen hergeleitet ist. 


(Received: October 10, 1960.) 


A Boundary-Layer in a Non-Newtonian Fluid 


By Joun R. JONES, Swansea, Great Britain?) 


1. Introduction 


The boundary-layer theory was developed by PRANpDTL [1]?) for fluids of 
constant viscosity, initially for the case of laminar flow, and later was extended 
to include compressible viscous fluids and turbulent flow in the boundary- 
layer. In recent years experiments have shown that several phenomena obser- 
vable in fluids are not predicted by the classical theory of viscous flow. This 
has led to the formulation of more complicated rheological equations of state, 
either deduced from a molecular picture or microscopic model of the material, 
or based on idealizations of simple experiments. The invariant forms of 
rheological equations of state for a homogeneous continuum, suitable for 
application to all conditions of motion and of stress, have been discussed ge- 
nerally by OLDROYD [2]; restrictions are imposed on their form by the physi- 
cal conditions that the rheological properties they describe are independent 
of the frame of reference and independent of the motion of the material 
element in space. With the growing interest in non-Newtonian behaviour it 
seems that one should now, if possible, extend the boundary-layer theory so 
as to include a wider range of fluids. In the present paper the effect of a 
variation of viscosity with rate of shear and of a normal-stress coefficient 
(representing ‘ cross-viscosity ’) are examined in a simple type of two-dimensional 
boundary-layer. For reasons of mathematical convenience the case which is _ 
worked out in detail, by way of illustration, is one in which the viscosity co- 
efficient is a linear function of the rate of shear, the normal-stress coefficient 
is arbitrary and the mainstream velocity is proportional to the one-third 
power of distance measured from a stagnation point. It is of interest to note 


') Department of Applied Mathematics, University College of Swansea. 
*) Numbers in brackets refer to References, page 343. 
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that some of the effects found here are qualitatively the same as those found 
when the compressibility of the fluid is taken into account (see, for example, 
_SCHLICHTING [3)). 
= 


uae 


2. The Equations Governing Non-Newtonian Flow 


__ There are various theories regarding different types of rheological behaviour. 
The simplest of all is based on the assumption that, apart from an isotropic 
pressure, the stress in an element of incompressible fluid at any instant is 
determined uniquely by the instantaneous rate of strain. In this case the 
‘theological equations of state are of the form 


‘ Pir = — Poin + Pix» (1) 
. here the partial stress tensor p;, is an algebraic function of the rate-of- 
; strain tensor hea 

i Ck =O us cy o) (2) 
“ 


at the same position and time only, and # is an undetermined isotropic pressure 
which varies with position and time in a way determined by the equations of 
“motion. Here v, is the velocity vector referred to Cartesian coordinates x; and 
-6;, is the Kronecker delta. On the assumption that the material is inherently 
‘isotropic and that the stresses can be expanded as power series in the rates-of- 
-strain, Rivrin [4] has shown that the most general possible equation of the 
_type (1) may be written in the form 


< 

f Dip = — POR +2 CAS Ci; Se» (3) *) 
_where the coefficients 7 and ¢ are arbitrary scalar functions of the flow in- 
variants 

* 1 

; I, =i) 1s = [(6:5)? — in Gcal I, = det (¢;,) 5 (4) 
£ 


“subject to conditions (such as 7 > 0) dictated by thermodynamic principles. 
‘We can consider equation (3) coupled with the incompressibility condition 
I, = 0 as a useful prototype set of equations of state for detailed study, on 
“the grounds of mathematical simplicity, without accepting RIVLIN’s con- 
tention that his assumption of a power-series expansion will be generally valid. 
A formula of the form (3) in which the coefficients 7 und ¢ are physical constants 
has been discussed by REINER [5]; REINER considers the material to be com- 
pressible and accordingly the quantity p is a function of the density i. e. of the 
dilatation. 

3) The usual convention of summation over the values x, y,z applies to all repeated Latin 
_ suffixes except %, y, 2. 
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Referred to a suitably chosen orthogonal Cartesian frame of reference — 
x,y,z, the components of velocity in any two-dimensional motion may be | 
taken as 


u=ulx,y,t), v=v%,y,t), w=O0. | (5) 
The equation of continuity | 
Ou Ov | 


is the condition for the existence of a stream function p(x, y, ¢) such that 


“=—-—, V= 


Ow Ow > 
iy ge | 


In terms of y the non-vanishing components of the rate-of-strain tensor eé;, 
and of ¢;; é;,, are 


eee O2y Be! Oy dy oe Ory 
ae “On Oy? “HU las ait yy “Ox dy’ 
Oy \2 1 (dp dy \2 
A aT ( Ox eal 4 ae ¥: a 


On substituting these values in (3) it is found that the stresses are 


foe 7 a , Pry mel (2 se] , 
by =—P* + 20a, b= — P80 [4 (oe) (8) 
(Se oe) |) t= Pa =O, 
where 
pare [2% + (Oe Sey, 9 


The shear stress #,, is independent of the normal-stress coefficient and is the © 
same as if the fluid possessed a single variable coefficient of viscosity y, and 
the normal stresses p, ,, £,,, P., are in general unequal. 
The stress equations of motion 
Dog. OL) = dye 
Q ( Dt : ie ) = th (10) 


Ox; OX, 


(where D/Dt = 0/0t + v; 0/0x,, @ is the density of the fluid, and U is the potential 
of the body forces per unit mass) require 


or so) +05 4 2 = om x ) 


oy Ox Ox 
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4 * 


pop (%) +0 +e clei borte oe lhearea Ch veserd ae 


where p* = p*(x, y, ¢). Eliminating p* from equations (11) and (12) we obtain, 
as the equation determining y, 


. 0 [2m bie eI Rees nee om] 


ee” fake 


; ot 0(x, y) ax Oy? On Oy? 
O*y 
= 2 il ie 5 (1 Ox Oy ), =) 
_ where 
. ney. Ot 02 
Vaare oye cays: (14) 


_ It is seen that equation (13) is independent of the normal-stress coefficient ¢, 
_ and so too are the boundary conditions (taken to be those expressing no slip 
at a rigid boundary) on y. Thus the velocity distribution for a general two- 
dimensional motion is the same as that in a purely viscous fluid characterized 
_ by asingle variable coefficient of viscosity 7. 
; The flow invariants in two-dimensional flow are 


L 20, R=—- 2; =0, (15) 
3 where 
: yo [+ (abh + (SE - ST 9 


the index 1/2 is intended as the positive root throughout. It is easily seen that y 
has a direct physical significance as the rate of shear at any point. It follows 
- from (15) that the viscosity coefficient and the normal-stress coefficient are in 
the most general case arbitrary functions of y. We shall in what follows confine 
- attention to the particular class of fluids characterized by an arbitrary normal- 
stress coefficient and a viscosity coefficient which is a linear function of y, Le. 


n= (l—ay), o=Cly), (17) 
_ where 7) and « are physical constants. The viscosity coefficient 7 is of necessity 
positive and this requires 

N20, L-aySod (allx,y,?#). (18) 
- Also, for physical interest, the shear stress must be a monotonic increasing 


function of the rate-of-strain, and we shall therefore have to restrict attention 


to regions in which 
No20, 1-2ay>0. (19) 


_ This restriction (which is stronger than (18)) imposes a limitation on the type 
of two-dimensional flow of these simple prototype fluids that we can study. 
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The equations (11), (12) and (13) reduce to | 
D oy Opt eS Uae 
e ol dy T Teg Ti8 oy 


4 na 2 ie (y ce iy ly (Se * ol}: ss 


D [op op* _—_—s- OV? 
alee ay. 18 Oe 


~ me fele (oe ~ Sl +2ap (ae )} 


Oviy O(y, Viv) ] _ 4 
[ Ot “ O(x, Ae =o Vi 


alte EIDE B +4ahs Hl). ea. 


where the potential function 9 U has now been incorporated in p*4). When | 


accelerations are zero or are negligible and the viscosity is invariable (« = 0) 
equation (22) reduces to the biharmonic equation. 

The rate at which work is being done per unit volume in changing the shape 
of an incompressible material is, in any motion, 


D= Pip ix, (23) 


independent of the rheological properties of the material. In any two-dimen- 
sional motion of an incompressible fluid characterized by (3), ® = 2 e;; e;,, 
unaffected by ¢, and when 7 is given by (17) we have 


D=2(1—ay) y?. (24) 


3. The Boundary-Layer Equations 


The arguments leading to a boundary-layer theory are very much the same 
as the familiar ones of PRaANpTL. It may happen that the typical length 
appropriate to the coordinates x and y are different, they may be written d and 
_e€d, with e < 1 without loss of generality. If V is a typical wu, we define dimen- 
sionless variables 


%=K'd, y=y'led, tV=t#'d, | 
p=pweVd, u=u'V, v=0'eV, (25) 
pr =p oV®, yed=y'V. | 


4) In what follows, therefore, — Pyx» — Pyy, — Pzz denote the excess normal pressures, in diffe- 
rent directions, above the hydrostatic pressure. 


RETO 
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_ These non-dimensional quantities are, with the possible exception of #’, all of 
_ order unity. In terms of the reduced variables equations (20) and (21) become 


af 1D) Ow’ Op’ 0 O2y’ O2y’ 
—2R _ +) oa bs (3% 2 Oy 
Som) eral Oy’ T oy Oy’ \ Oy” aya 
e ie Ga 0 O2y’ Oy’ 
1 2 y , y ee) wy 
<8 Q{2 ex’ (v Ox’ Ov’ xr Oy’ ly S. o alle (26) 


D [ey’\ op’ 0 (0% Oey’ 
2 R 9 (57) |=#4r( y 2 ¥) 
3 oe pete Tay? | == pa apt t Po gx 


+ 20ffy (Sm —e2m)] anti s2%|}. cn 


; where R, and Q are the dimensionless parameters defined by 


; _eVvd av 

Z Boe ch ins aa (28) 
ro (( Ow. _ 2 Oy’ \? 2 (_ oy’ \2712 

, Y= [(Ge— aan) +44 (arar) | (29) 

and 

| a 300 ao ue, 

; De oF oe He: (30) 

. 


The orders of magnitude of the groups of terms in each of the equations (26) 
_ and (27) are (from left to right) e? Ro, 1, e~* Q and e« Q. Interesting cases arise 
_ when two or three of these four quantities are of the same order of magnitude 
in part of the region of flow, while the remaining quantities are smaller. If 
' ¢? R, is approximately unity and if both e~1Q and e Q are small, we have a 

Prandtl boundary-layer for a fluid of constant viscosity 7. In the fluids now 
under consideration the possibility arises of a region in which e? Ry and e~* Q 
are approximately unity and ¢ Q is small; this may happen when R, is large 
and Q is small so that Q Ri? is of order unity. In such a region of flow (a 
- boundary-layer) certain terms may be neglected in the equations of motion 
~ (20) and (21). We shall be concerned only with forward flow in the boundary- 


layer and therefore may take 0?p/0y? < 0 so that y’ = — O*y’' /Oy'*, approxi- 
mately. Neglecting terms of order ¢*, equations (20) and (21) become 

pee or ee, OW 8 (Se A oP Lo 31 
a dhe og 0x ?0 Oya Yet (a) : o oy ; (31) 


_ where % = Mole. The normal stresses p,,, Py, and p,, can be written 
es > (Ov \? 
Pes =Pyy=— 2» bar =— 0° —0(55) - (32) 


The stress component #,, is here a function of x, y and ¢, whereas for Newtonian 
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fluids p,, is a function of x and ¢ only. The shear stress in the boundary-layer is 


r 


2. o2 
t= piy=— 1 (1+ 05%) 5% (33) 


| 


and the rate of dissipation of energy per unit volume in the boundary-layer is _ 


O?y\ ( O?p\2 
@ = o(1 + 058) (G58) - am 
The inequality (19) reduces to 
gs 1+2a5%>0 (all x, , 2). (35) 


In most parts of the fluid the viscous stresses are negligible, and just outside 
the boundary-layer the flow is that of an inviscid fluid defined by a mainstream 


velocity u(x, #) and a pressure (excess above the hydrostatic) #,(x, 4). The | 
continuity of #,, at the outer edge of the boundary-layer requires that — 


p*(x, t) = p(x, t), and, in the case of steady motion, the equation determining 
wy becomes 


op dp oy Oy du, dy F) Gar (36) 


dy Ox dy Ox dy? “1 dx ~ "0 Qys — "0% By \ Oye 


The boundary conditions under which equation (36) is to be integrated are | 


the same as in PRANDTL’s theory. 


A search for solutions of equation (36) in which the velocity profiles are — 


similar in different sections reveals that the only possible forms of u,(x) and 
(x) (the boundary-layer thickness) under these conditions are 


fh, OCF OOS AMS. (37) 


where x measures distance along a rigid boundary from a stagnation point). 
Taking the mainstream velocity as k x1/3 (k > Q) it is convenient to write 


gibad sae 2s F(é) Phe. 2k \ x18 y ’ (38) 


3% 


so that equation (36) reduces to 
F"(6) + FQ) F"@) <K F"® P"® = 4 (FE) — 1, (39) 


where a prime denotes differentiation of a function of a single independent 
variable and the dimensionless parameter K (which is proportional to the 
dimensionless quantity QR"? and is, therefore, of order unity in the boundary- 


5) Cf. GoLpsTEIN’s [6] method of finding similar solutions of the Newtonian boundary-layer 


equations. When Q = 0 an infinite number of similar solutions exists, defined by u, 0c x™, 
6 co ~~ 1/2(m—1), 


poe” 


Vol. XII, 1961 A Boundary-Layer in a Non-Newtonian-Fluid 335 


- 


Bayer) is defined by 


8 k® \1/2 
; K=(3,) (40) 
Equation (39) is now to be integrated with the boundary conditions 
‘ Oven F(Ojm 0. fiio0) = 1. (41) 


We shall work out the details of this particular case in §§ 4, 5. 


4 4. Approximate Solutions Based on the Momentum Integral 


F Integrating equation (36) across the boundary-layer we arrive at the 
momentum equation, namely 


5 5 
ae? f. Op \ Oy duy oy 
ig [ (at ales dy ae f (m+ $e) 

0 0 
4 oT et) ane 
7 A Roa cdets eal alent rete (42) 
4 If, as in PRANDTL’s theory, we define a displacement thickness 6,(x) and a 
- momentum thickness 6,(%) by 


4 : i 
Be anfo-s)o. a-f0-Do. 
3 0 9 


and express (42) in terms of 5,, 5, and T, (the shear stress at the wall), we are 
able to reduce the momentum equation to the same form as for a Newtonian 


fluid: 


d do 
ty Ger (1 + 285) + wig = (44) 


; We now take u to be approximately a polynomial of the fourth degree in the 
- dimensionless variable X = y/6 whose coefficients are functions of x only: 


; u = uy(f, X + fg X? + fy X* + fy X%) « (45) 


Using the same boundary conditions as in POHLHAUSEN’S [7] approximate 
treatment in the case of Newtonian flow we find that 


1 1 
h-2+42, h=-Za4, f=-2+q4, h=1- Gi (46) 
where roe 
au 4 au Cie 
EM ya _ (1 48%), 4 % Bao, (47) 


6) It is appropriate to take Ry = k x4!3/p, and Q = « k/x?!3 in this problem. Near the stagnation 
point (* = 0) the boundary-layer theory ceases to be applicable (as in the case of PRANDTL) since 
the assumptions that R, is large and Q is small are not satisfied. 


1 
~ | 
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(The equations familiar in POHLHAUSEN’S treatment of the ordinary boundary- 
layer are obtained by putting « = 0.) Hence we obtain for the velocity profile — 


=a BCO ea CN: (48) 


where 


H(X)=2X—-2X84X*, 6(X)=+xX(— x). (49) 


The forms of H(X) and G(X) are exactly as for Newtonian fluids and have been — 
exhibited graphically (see, for example, GOLDSTEIN [8]), but here A depends © 
on both the dimensionless parameters Ry and Q. . 

The boundary-layer thickness must be of the form 


i= 8(a3) =. (50), 


where # is a constant; the parameter / then has a constant value throughout | 
the boundary-layer and equation (47) reduces to 
3 f?-—6AB+ KA(A4+12)=0. (51) 


On the basis of this approximate analysis, the displacement thickness and the — 
momentum thickness are given by 


See, (S2)" a8, 6,= Be ( at) xls (52) | 
where 
She kt 2 frail I 5 
fay BOD, RR (8) Se ee) 


equations already familiar (with a different meaning of 4) for Newtonian flow. _ 
The values of f, and £, are easily found when 4 and f are known. 
Equation (44) now reduces to 


(548 — ios Ors 72) p> — 280 (a Ae+ kuz y= (54) 
3 18 3 4 

and equations (51) and (54) may be solved numerically for a given value of 
the parameter K. The values of A, B, 8, and f, for various values of K are | 
included in Table I. In order to avoid ambiguity in the values of A (i. e. of ) 
we must consider only the range 1:83 < 4 < 12, and then the results derived 
by the above approximate method are unique. It follows that we must cer- 
tainly have K < 0-41. The restriction (35) implies 


K 
il  —- 
B 


and is consistent with the range of values of K < 0-41. 


(2+ 4) >0 


SO 


.* 
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Table I 


(3 8/8 r9)}!2 a2! Cy 


0 05 10 -14 -12 -08 -0-4 
Oy 


Figure 1 Figure 2 
Showing how the coefficient of friction varies 
with the parameter K. 


Showing how the velocity profile across the 
boundary-layer varies with the parameter 1G 
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The velocity profile across the boundary-layer corresponding to each of the) 
values K = 0-4, K = 0 and K = —1-4 is shown in Figure 1. Comparing the 
columns of Table I it is seen that the displacement thickness (proportional to| 
f,) and the momentum thickness (proportional to f,) both decrease as K in-| 
creases over the whole range of values of K considered. (The values of 8, which 
has a less direct physical significance, do not reflect the trend in boundary- | 
layer thickness so well.) 


The skin-friction coefficient C; is given by 


Gora ap te )ikoay Coe laa) ears 50), 


and the values of (3 k3/8 v9)1/? x?/8 C, for corresponding values of K are included — 
in Table I and are exhibited graphically in Figure 2. It is seen that as K 
decreases from 0-4 to — 1-4 the quantity (3 k/8 y9)"? «2/8 C, increases almost 
linearly with K. 

The normal stress p,, is given by 


Peet Pr= — sy nae (A + 12) - 62 X+9(A—4) X46) AHP, (57) 


where X varies from 0 to 1 across the boundary-layer. At the outer edge of the 
boundary-layer (X = 1) we have 


(p22 + Pi)y 9 =0- (58) 
At the wall (X = 0) we have 


(Oh tae = ee (A+ 12)2. (59) 


Therefore normal stresses perpendicular to the plane of flow are required in 
the boundary-layer region near a rigid boundary in order to maintain a two- 
dimensional motion, the sign of these stresses depending on the sign of ¢. The 
magnitude of (p.. + ~1),.9 depends on the parameter K; (22+ Prjy=o 54 
creases with K and the value at K = 0-4 is roughly 2:4 times that at K = — 1-4. 


5. A Method of Evaluating the Parameter F’’ (0) 


The accuracy of the modified Pohlhausen method of § 4 can be checked by 
comparing the results with those obtained by a more direct method of inte- 
grating the non-Newtonian boundary-layer equations. A method recently used 
by MEksyn [9] in solving the boundary-layer equations for Newtonian flow 
can be extended for this purpose. It is assumed, in the range of values of the | 
parameter K under consideration, that the boundary conditions (41) determine — 
F(é) without ambiguity and that the value of F” (0) is a continuous function of 
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K. The inequality (35) may be written as 
4 1—K F"(0) >0 (60) 
and there appears to be no other restriction on the magnitude of K. 
__ It may be shown, by inspection of (39) and (41), that we must seek a solution 
of the form 

PONE) AMe) exe ey; (61) 


where the positive function N(€) is given by 


€ 
N@ = | Fé) a, (62) 


and M(é) is a slowly varying function; a solution of this type, satisfying the 
boundary conditions at € = 0, ensures the correct form of the function at 
& = co. We observe that the function F'"(€), as defined by (61), is of precisely 
the same form as for Newtonian flow (K = 0), except that here the functions 
M(é) and N(é) depend on the parameter K. Hence the steps of MEKSYN’S 
analysis may be followed closely in the non-Newtonian case also, and the 
‘reader is referred to MEKSYN [9] for fuller explanations. In the notation of 
‘MEKSYN we write for the functions F(£), M(é) and N() the power series 


Fe-e La MO-DAe: MOP Lae, (63) 


with constant coefficients @,, 6, and ¢,. 
_ We require to find the value of the parameter F”(0) that permits the 
boundary conditions at infinity to be satisfied; and we can use the condition 


F'(co) = ‘: M(é) e-N@ dé =1. (64) 
; 0 
To evaluate (64) we put 
q o = N(é), (65) 
which implies 
c=) Ao (66) 
n=0 
where 
il qr-1 s 
Ay = oy [eer (Oo + aE + Ob +o) Boma: (67) 


Flt may be shown that (64) reduces (in terms of gamma functions) to 


co 


lim dant (252) wttat, (68) 


G10 
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where | 
yi ih qr 2 j 
0, = xi {Gor (Po + E+ Dg +) | 
x (Cn si te é Ss oe & a shaogo . (69)} 


0 


The coefficients A, and d,, expressed in terms of the parameters F”(0) | 
K, are 


6 \1/3 (an Ds 1/3 
alee) On A, = (2) A A, = == (4) , 
a~* L2 
= Ie 
a) LL? | a@\2/3 
= 16 L? — 684 L + 643) , 
ie 5760 is) (2 eon) 
3 —272 77 a-8 L4 lla->*KLS5 
A ae a 1/ aa a za % (70) 
5 6 140 13824 960 
3 a-4 K2 6 se ee 
* 200 - 2240 g 
¥ a-* [3 a~2 K L4 9a-7 Ls 3a-*K L8 
cmte Seg tl 585; 8192 1024 
63 a-> K2 17 9 a-* K3 L8 3 a-3 K4 L$ 
12800 + 1400 512 : 
i 5 6 \2/3 
dg = (6 a dy eee (2 )™, 
3a-2L2 
dg= To (2 L-17), 
aes 6\1/3 (a 91 a-3 L3 Das kK Ee a) 
a (=) (3 5184 144 _ 60 
Pike: “sell aL EPs 935 a-4 LA 43 a-3 K LS 
eels (= 36 126 248832 3456 
93a-2 K2 Ls ee) 
~ 2880 " 366 ; 


*) Meksyn’s results for the constants d,, dy, d, and dg, in the special case when two of his para-_ 
meters («% and 4’) vanish, can be obtained Been ie above expressions by putting K = Oe ee) 
The corresponding yesulls for the constants dp, d, and d,, obtained by putting K = 0, L = 1 in 


(71), are 
1 6 \2/8 / 25 935 a—4 
dy = — (6a%))'2, a--(- if 
a 252 248832 


4 -(2)" (4 169 a~8 43225 a7 
‘ 6 20736 35831808 } ? 


- 
=> 
z 
_ 
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eel a S18) | 297 Ligh Boas KGL =) 2078 L! 
“iq 160 1120 140 > AOS6 
5Sla-> KILLS 9a-4K2L?  3267a-2 K3L8 81a-2 K4L? 

: 2560 = 160 22400 & 256 : 
7 ie 4 oles OSs aeL* 2080" %R 15 

a 6 25920 34560 1152 (71)?) 

is 7a-1 K2L5 43225 a-7 L8 21 a-! K2Lé 
480 35831808 128 
. he 5369a-* K L? 13741 a—5 K2 L8 8591 a—4 K3 L 
; 995328 829440 172800 
3629 a-3 K4 L100 245 a-? KS [1 
30720 1024 : 
where 
g= 0) “and “b= (i a Kj+1. (72) 


__ The following procedure is now adopted to complete the calculations. For 
a selected value of the dimensionless parameter K we assume a trial value for 
the parameter a. The series in (68) is, in all the cases considered, divergent as 
~ tends to unity, and we must apply EuLEr’s transformation (see, for example, 
Harpy [10]) repeatedly until a convergent expression is obtained. An adjust- 
ment in the value of a can then be made, and the process repeated until (68) is 
satisfied. The results derived by the extended Pohlhausen method of § 4 are 
used to give a first approximation to a, and the calculations are continued 
with steps of a = 0-005 until the result F’(oo) = 1-00 is obtained. The numerical 
calculations are heavy and, for the sake of shortness, only the final results are 
summarized here. 
The constants f, and f, are 


f= f0-F@le, m= fF@U-F@Me. — (73) 
Hence : ‘ 
B= f EMU) ec dé, (74) 


which follows on integration of the first equation of (73) by parts with the 
assumption of the property . 
Hee Piel) ee (e) 0 (75) 


which differ from the expressions in MEKsyn’s equations (9.7) with the particular values « = 0, 
WV =0,A = — 1/2. A simple misprint accounts for the discrepancy in the case of dy, and the dis- 
crepancies in d, and d, are of no more than about 0-1 %, so that Mrxsyn’s final results 
are not in question. 
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Applying the transformation (65) we obtain 
Soe = WPL sans 7 6 

B= lim Se, T"( 3 )x ’ ( y 


where 
é, = A, d, + Agd, 1 ++: + Anyi dy, (77), 


and to sum the expansion for x = 1 we use EULER’s transformation. 
To find the constant fh, use is made of the relation 


[1-K F"0)| FO) 


-; f/0-FO@le+s5 /FOU-FOlae=-LA+Fhs (78) 
0 0 


which follows on integration of equation (39) with respect to & with the as-- 
sumption of the property (75). 
The final results are exhibited in Table IT. 


Table II 


(3 k3/8 Yo)i!2 42t3 Cr 


A | es 


0-4 0-74 0-31 0-85 
0 0-80 0-35 0-93 
—1:-0 0-90 0-42 1-08 


6. Conclusion 


The values of B,, By, (3 k?/8 v9)1!? x28 C,.as found by the two above methods: 
of solution agree fairly well and one is therefore led to regard with confidence 
the above results based on the adaptation of the momentum equation. 

Although the quantitative results in §§ 4, 5 have been obtained for a parti- 
cular type of non-Newtonian fluid (i.e. for a particular form of the viscosity 7) 
and for particular external conditions, we may expect them to be qualitatively 
valid much more generally. The variation of viscosity with rate of shear must 
be expected to decrease or increase the parameters 6,, 6, and C y according as the 
viscosity decreases or increases with the rate of shear, and any normal-stress 
coefficient will not be expected to affect the velocity field. | 

I should like to thank Professor J. G. OLDRoypD for most helpful discussions 
and suggestions throughout the progress of this work. 
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Résumé 


Dans ce mémoire, on examine l’effet d’une variation de viscosité avec la vitesse 
de déformation dans une couche limite a deux dimensions. Pour des raisons mathé- 
‘matiques traitables, l’exemple que l’on étudie en détail est celui dans lequel le coéf- 
ficient de viscosité est une fonction linéaire de la vitesse de déformation et la 
-vélocité du courant principal est en proportion du tiers de la puissance de distance 
‘mesurée d’un point de stagnation. On peut noter que quelques-uns des résultats 
“qui se trouvent ici sont qualitativement les mémes que ceux qui se trouvent quand 
‘on tient compte de la compressibilité du fluide. 


" (Received: October 24, 1960.) 


Energy Flows ina Vortex Tube 


By ALan J. ReyNnotps, Cambridge, Great Britain’) 


In a recent paper [1]2) the present writer considered the dynamics of the 
‘turbulent flow of a compressible fluid using an order-of-magnitude analysis. It 
‘is now proposed to extend this analysis to include energy fluxes in such vortex 
flows and, in particular, to give an explanation of the energy separation which 
is termed the Ranque effect. 


1. The Ranque-Hilsch Vortex Tube 


It will be remembered that the Ranque-Hilsch vortex tube is a device which 
is capable, although without moving parts, of dividing a homogeneous inflow 
of gas into a cooled stream and a heated stream (reference being in each case 


1) Cavendish Laboratory. 
2) Numbers in brackets refer to References, page 356. 
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to total temperature). The usual form of the device is as follows. The chief 
component is a tube of constant circular cross-section, its length several times 
its diameter. Near one end the tube wall is pierced to allow the entry of fluid! 
(usually air), the inlet ports being so disposed that a swirling flow is set up) 
within the tube. The fluid escapes at the ends of the vortex chamber, some} 
through an orifice plate near the inlet, the rest through a second throttling; 
device at the other end. Under certain operating conditions (notably, for’ 
reasonably high supply pressures and flow speeds) it is observed that the air’ 
emerging through the orifice near the inlet has a lowered total temperature and_| 
that the air leaving at the other exit has a higher temperature than the air: 
supplied; temperature differences equivalent to several times the maximum. 
dynamic head can be obtained. : 

The device and the effect it produces are associated with the names of two) 
men: RANQUE, the French engineer who first fashioned it in 1931, and Hirscu, . 
the German physicist who first studied its performance in detail and applied . 
it to the liquefaction of gases in 1944. Very many other investigations of the 
device have been made since their early work. WESTLEY [2] has given a 
bibliography of work before 1952. Since then further performance data have 
been published by OTTEN, SPRENGER, and WESTLEY; results of internal flow 
measurements by HARTNETT and ECKERT, Lay, SCHELLER and Brown; and 
theoretical treatments by DEISsLERand PERLMUTTER, HIwn2E, Lay, PENGELLEY, 
and Rott (([3] to [12]). Several other studies have been made in this period, 
but they are less important or less easily accessible. 


2. Qualitative Theories of the Ranque Effect 


Early discussions of the energy separation within the vortex tube considered 
inviscid, non-conducting fluids: the first theory put forward by RANQUE was 
of this type. However, it soon became apparent that transport properties must 
be considered. A later theory of RANQUE, supported also by Hiscu, was that 
the outwards migration of energy was due to radial work transfer. They pictured. 
a ‘free’ vortex being converted by viscous shearing forces into a ‘forced’ 
vortex and a consequent outwards transfer of kinetic energy. FULTON [13] 
carried forward these ideas of vortex conversion by noting that inwards 
conduction of heat would counteract the outwards work flux. He found that 
a net outwards flux occurred in laminar flow but estimated that it was not of 
sufficient magnitude to produce the observed energy separation. 

It thus appeared that turbulent mixing was a more likely agent of energy 
transfer. This conclusion still seems valid today: it is supported by the recent 
calculations of DEISSLER and PERLMUTTER [9]. | 

At first the role of turbulence was thought to be qualitatively the same as 
that of molecular transport, although giving rise to very much greater energy 
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fluxes. The first discussion in which a contrary view was taken is that of 
_KNOERNSCHILD [14]. He pointed out that a compressible fluid lump moving 
rapidly through a pressure gradient must experience temperature variations 
on expanding or contracting. Thus the temperature distribution will move 
_towards that related to the pressure gradient by an adiabatic law. At present, 
-KNOERNSCHILD’s concept is the most widely favoured explanation of the 
‘Ranque effect, although a work flux contribution is also thought to be impor- 
tant by some workers. 
The ideas just discussed form the main stream of thought concerning the 
-Ranque effect. But several suggestions have been made which do not fit into 
this evolutionary pattern. The most challenging is that put forward by SPREN- 
GER [4], but seemingly due originally to ACKERET, that organized unsteadiness 
(as distinct from turbulence) might produce the energy separation. This hypo- 
thesis is based on analogy with resonance tube phenomena. The present writer 
_has already put forward his reasons for feeling that this mechanism offers a 
less satisfactory explanation of the Ranque effect [15]. 
| Summarizing: it seems probable that the energy separation of the vortex 
tube is the result of turbulent activity within the vortex chamber. Our task 
is to discover what are the mechanisms of turbulence by which the energy 
_ fluxes are set up. 


3. The Energy Equation 


| If we neglect molecular transport properties and body forces, the energy 
equation can be written 


0 0 
L(g h) +V- (Qhou) = <P, 


where hy = h + u?/2 is the total enthalpy of the fluid, the sum of the enthalpy 
and the kinetic energy per unit mass. Taking mean values with respect to time 
and requiring that the mean motion be steady, we have 


V- (of u) =0 
which gives 


1) 0 
py 1 EH ho) Ai 9; (2 w ho) =0 


for axisymmetric mean flow. We have 


ouhy=(R+o') (U +4’) (Ho + ho) 
and ERT 
Hy =H +5 (U2+V?+ W2+ 42+ v2 + w’), 


Kg=h'+wU0+uvV+u0'W. 


346 Aan J. REYNOLDS AME 


4. Order-of-Magnitude Analysis 


3 

We shall specialize the energy equation given above by introducing a 
system of approximation suitable for vortex tube flows. The scheme to be 
used is that set forth in Section 3 of the previous paper on the dynamics of 
vortex flows: 


As before, we specify the turbulence by 
uw, 0, w® ~O(B), 9’ ~ O(n), 


adding here h’ ~ O(a). 

We choose units such that 7, V ~ O(1), knowing from the previous ana- 
lysis that R ~ O(1) for consistency. We take also 4H, ~ O(1), knowing that 
energy separations are typically a few times the dynamic head. This assumption 
will be seen to be consistent with the others a posteriori. It will be assumed, 
of course, that a, B, 06, €, 7 are all small compared to unity so that terms of 
higher orders in them can be neglected. 

Introducing these approximations into the appropriate form of the con- 
tinuity equation, we find that 

é€ 
Wend: CB 
at most, so that 
yg CRU) + FARW) + 2S (rel u!) = 0618 67) 
O(e) O(e) O(B*!* ») 


is the approximate continuity relation for this case. 
Similarly, the approximate energy equation is found to be 


ih fy 0 1 <0 BORN: 
Peep iop ee) aR Wie) = oy " Hoe u’) 
O(e) O(e) O(B*!? n) 
Lo eee EO ert AK! ae 
a NR Wt ae a ne v) Pr Rte w’) 
O(a B*!*) O(B) O(Be/d) 


= O(B 6, Be, « Bt? 6). 


Only the leading error terms have been indicated. : 
This result relates the several important energy fluxes into an elemental 
volume. It can be cast into another form without altering the orders of the 


as 
2 
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_ errors if we use the continuity equation given above. Thus 


> Diy. s. if oe) eran je OTL 
a4 R OE olla ate 5 RU R) ou i 
20 TF an 18 Per 
sire rs AL ee a w’) 


where D/Dt indicates the rate-of-change following the mean motion. Here we 
_have a relation among the radial energy fluxes and the energy input to a lump 
_ of fluid moving through the system. 


5. The Important Energy Fluxes 


The simplified forms of the energy equation suggest that the significant 
energy movements within the vortex are due to four energy fluxes: 
| a) The heat flux Rw’ h’ is produced by turbulent mixing through the radial 
_ temperature and pressure gradients. This is the contribution whose full signi- 
_ ficance was pointed out by KNOERNSCHILD. 
b) The flux of total energy 0’ w’ Hy is associated with the ‘Archimedean’ 
correlation which was discussed briefly in connection with vortex dynamics. 
c) The work fluxes R V uw’ v' and R W w’ w’ associated with the two impor- 
tant Reynolds stresses, those acting on cylindrical surfaces. 
A rough estimate of the relative magnitudes of these contributions can be 
_ made using data from extensive tests on vortex tubes. These data were used in 
a similar manner in the previous analysis of flow dynamics [1]. The vortex 
tube within which the flow was investigated was about 7-5 cm in diameter and 
about 120 cm long. Total and static pressure and temperature traverses were 
made at several stations down the tube with various blockages at the ends. 

In [1] it is found that the results suggest that we may represent a typical 
internal flow by taking 


oe 


Pena Wri e: fy ke, 
50 2 200 ’ 100 ’ 100 


Thus we estimate the terms of the energy equation obtained above as 


eG) 0 LENG, pone, 
<5, (1 RU Hy) + (RW Hh) + — 4 -(r Hoe u’) 
0(1/200) 0(1/200) (1/1000) 
Lee ah 1.0 ea inc mre 
sal Yo lie od lsc gies Me ao ay T RW w’) 


O(a/10) 0(1/100) 0(1/400) 


1 (4 
Ol-s00 300) 
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We note that the Archimedean effect seems to be somewhat less important | 
that the other contributions. The work fluxes appear to be of the same order . 
as the heating/cooling effects (the first two terms). There is no reason to 
suppose that the energy flux due to the Knoernschild effect is more important 
than that produced by the radial work fluxes. The quantity « will have to be 
~ 1/20 if the thermal convection term is to equal the work flux in importance. — 

We next consider the physical processes corresponding to these important 
terms of the energy equation. In this way it will be possible to obtain a better _ 
idea of the relative importance of the several contributions to the radial energy © 
flow. We adopt now the so-called Lagrangian point of view and attempt to 
build up idealized models of the internal motions of the turbulent fluid in | 
order to see how the energy flows are influenced by mean-value gradients in — 
the vortex. 


6. The Knoernschild Effect 


KNOERNSCHILD observed [14] that it is not correct to think of the turbulent 
mixing of a compressible fluid across a pressure gradient in terms of the 
simple model used for an incompressible fluid — a fluid lump shifted unchanged — 
to a new position and there mixing with the ambient fluid. A lump of com- 
pressible fluid will not be transported through a pressure gradient with all its — 
properties unchanged and without interactions with the fluid through which — 
it moves. A work interaction with the surroundings will occur during such a | 
movement, the particle being compressed on moving radially outwards to a 
region of higher pressure or expanding on moving inwards. Thus the work 
interaction with the matrix of fluid through which mixing takes place results 
in the cooling of the inwards-moving particles and heating of those moving | 
outwards. 

It has generally been assumed that the convecting particles do not interact 
with their surroundings except by expanding and contracting. That is, it is 
held that they are very nearly adiabatic while in motion. DEISsLER and PERL- 
MUTTER [9] mention work supporting this idea. In any case, the utility of the 
theory does not depend on strictly adiabatic lump movements. The most _ 


recent derivation of the expression for w’ h’ in terms of mean-flow gradients, 
that of Hinze [10], takes account of non-adiabatic lump motions’). 
Assuming the mixing fluid to be a perfect gas, H1nzeE gets a result equi- 
valent to 
uh! = 


or y P or 


fee jeieH ae 


*) It should be noted that we are able in this way to take account of the effect of molecular | 
transport properties on the small-scale motions of the turbulence, even though viscosity and | 
thermal conductivity were neglected in the mean-value equations with which the analysis was 
begun. 


tp 
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For this case we have simply H = c, T, where T is the absolute temperature. 
Here ey is a positive coefficient expressing the activity of the turbulence — an 
_eddy-diffusion coefficient. This result shows clearly the tendency of the radial 
distribution to move towards that defined by an adiabatic law applied to the 
pressure dustribution. 


Using the radial equilibrium condition: 
oP y? 
es ob Ol) 


- 


an alternative form can be given. 


rh! OH V2 
ks in tu (5 - +): 


The simple dependence of the radial flux on the local centripetal acceleration 
is thus made clear. This result is valid only for a perfect gas, even though the 
nature of the fluid does not enter explicitly. 

. When the radial temperature gradient is less steep than that defined by the 
adiabatic law (as was the case in the vortex tube flows investigated) the heat 
flux due to this effect is everywhere outwards. The influence of this term may 
_be expected to be strongest in the core of the vortex where centripetal acceler- 
_ations are highest. 


7. Archimedean Effects 


So far the influence of the heat flux, R w’ h’, on the overall energy transfer 
has appeared straightforward. But YUDINE, in the discussion of HINZE’s paper 
on energy transfer in the atmosphere [10], pointed out that his treatment is 
not quite complete. In addition to the vertical transfer due to non-selective 
mixing by turbulent stirring (that considered by Hinze), there exists in the 
atmosphere an upwards heat flux set up by Archimedean forces. A counter- 
part of this phenomenon must be expected in vortex flows. A buoyancy effect 

: (0' w’) did emerge in the simplified energy equation for vortex flow. But this 
term does not wholly account for the effects of Archimedean forces. This can 
be seen as follows. It is apparent that a particle lighter than its surroundings 
will tend inwards under the action of buoyancy forces. But a deficiency in 
density can be thought of as resulting either from a lower than ambient 
| pressure or from a higher than ambient temperature. The first of these would 
seem to contribute primarily to the energy flux associated with the Archimedean 


correlation function, 0’ u’: its discussion will be deferred. The second possi- 
bility, low density corresponding to a temperature above ambient, is of 
immediate interest, for it contributes directly to the correlation function w’ h’ 


(as well as to 9’ w’). 


| 

4 
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not hard to see. Warmer lumps will tend inwards so that an additional heat flux | 
is set up, one which exists even if the mean-value temperature distribution is © 
adiabatically related to the pressure distribution and which is inwards in all — 
circumstances. The estimation of this buoyancy contribution in comparison — 
with that suggested by KNOERNSCHILD and studied by Hinze and others is 

not easy. We shall not study this difficult problem in detail here: PRIESTLEY © 
and SWINBANK [16] have discussed the analogous atmospheric case. But it is | 
impossible to leave this matter without reaching some decision, however 

arbitrary. We shall then assume that this Archimedean heat transfer is masked 

by the Knoernschild effect so that the heat flux associated with the correlation — 


| 
The general effect of this contribution to the heat flux in a vortex is . 
i 


function w’ h’ is outwards almost everywhere in the vortex. Support for this 
assumption is given by the following considerations. First, the Archimedean 
contribution to the heat flux can be of decisive importance only when the | 
temperature distribution is near that prescribed by the adiabatic law applied 
to the pressure distribution. And, in the vortex tube flows studied, this was — 
seldom the case. Second, it seems evident on considering the temperature 
distributions found in the Ranque-Hilsch tube that the outwards fluxes near 
the core do, in fact, overpower this additional, opposed heat flow. Lastly, the 
effect is of secondary importance in the analogous atmospheric situation. 

We have tried to separate out three influences of compressibility on energy © 
transfer in a vortex flow — the Knoernschild and two Archimedean effects. — 
The separation is a somewhat arbitrary one: each of the effects must in reality 
interact with the others. This approach does display some of the possible 
phenomena even if it does not clearly define the relative importance of the 
contributions. 

We now turn to the buoyancy effect which did emerge explicitly from the 
energy equation. As this contribution has not been dealt with previously, it 
will be given fairly detailed consideration here. Later an attempt will be made 
to relate the Archimedean correlation function, o' wu’, , to the mean-flow gra- 
dients. For simplicity, we restrict consideration at that point to density per- 
turbations induced by non-uniform pressures in the turbulent field. The 
methods used can be applied to other models. 

As was pointed out in the discussion of flow dynamics, the correlation 


function 9’ u’ may be expected to be positive in a vortex in a compressible 
fluid: lighter lumps of fluid will tend inwards under the action of buoyancy © 
forces. Then the associated energy flux, ow Hy, will be everywhere radially 
outwards. Surprisingly, the overall effect of this outwards flux, as indicated 


by the term ni 0H,/Or, seems to be to cool the fluid throughout the tube _ 


(since, except very near the tube wall, 0H,/dr > 0 in almost all of the situations 
investigated.) 
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_ This conundrum can be resolved as follows: the overall effect of the energy 
flux is to redistribute energy in the various regions of the tube; but, when 
continuity requirements are taken into account, this results in a reduction in 
‘the energy per unit mass of the whole of the moving fluid. Consider a simple 
example: purely radial flow. Continuity is 

au eu) =0; 


‘whence, for zero net source flow, we get 
z RU=-o'w'. 


We shall take o’ u' to be positive. Then the requirement of continuity is an 
inwards mean velocity, U=— 9’ u'|R. The energy equation derived from 
eeinuity is Fenwet Se az 


fame 
RU or Ry eer ny 


‘Thus we see than the rate-of-change following the mean motion is, in fact, the 
‘rate-of-change as the fluid moves inwards in the direction opposite to the 
energy flux. 

We shall next attempt to relate the correlation function g’ w’ to the mean 


flow as has already been done for the function w’ h’. Let us clarify the situation 
by relating the turbulent velocity fluctuation at a point to the perturbation 
velocity of a lump at that point. We write u’ = u, + w, where w’ is the deviation 
from the time-mean value at a point; w{ is the velocity of the coincident 
particle due to the Archimedean buoyancy force per unit volume, 0’ V?/r; and 
uw, is the part of the local motion induced by the surrounding turbulence. The 
term u; is random with respect to the fixed point but not so far as the passing 


eddy is concerned. Then 9 w’ = 9’ Uy. 
Our next task is to relate the velocity w/ to the buoyancy force producing it. 
It seems plausible to compare the inwards-moving lighter lumps to microscopic 
particles whose net drift through a gas is superposed upon Brownian motion. 
In the present case, the turbulence of the fluid supplies the random disturbance 
“to the drift due to the Archimedean force. In his theory of Brownian motion 
EINstEIN [17] has assumed that the drift speed is linearly related to the force 
on a particle, in fact, that the drift motion obeys STOKEs’ law. The results of 
his theory have been verified experimentally. Similarly in electricity: a linear 
resistance law — OuM’s law — is found to apply when the electrons’ drift speed 
is small compared to the speed characterizing the random thermal motion. 
By analogy with SToKEs’ law we relate a typical density fluctuation and a 
typical drift speed by 


, 


2 , 
0 VE oc Rey 


nt | 
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where FR é,, replaces the molecular viscosity, é4 being the eddy viscosity. The 
constant of proportionality is dependent on the scale of the turbulence. 
Introducing once again (as in [1]) the hypothesis that density fluctuations 

at a point are the result of the convection past that point of eddies in which | 
density gradients are maintained by pressure gradients holding the eddies 
together, we obtain 


"oo Ru’; Ay ¢ 
pocku p 


where K is the compressibility of the fluid. Then 


2 
(i 2 J d 
QUE, RK ean 


where é4 is a positive coefficient dependent on the kinematic aspects of the 
turbulence. 
Alternatively, introducing the radial equilibrium condition, we get 


eG) oP 
OU ok) on 


for the Archimedean correlation function. 


8. Work Flux due to Tangential Stresses 


To investigate the work fluxes we need not consider the details of the 
turbulent motion. The previous investigation of the dynamics of vortex tube | 
flows provides the means for the study of the work transfers. It was found — 
that tangential momentum conservation could be represented with reasonable 
accuracy by 

ess " | RW 4 Riri. 


v2 Or 


In its present form this result does not easily yield information about the 


stress Ru’ v' which gives rise to the work flux of interest here. An obvious 

simplification is the restriction to swirl distributions independent of axial 
position. Fortunately, this idealization has close counterparts in the vortex 

flows studied experimentally; and in particular it applies to the configurations 

for which the Ranque effect was especially large. The conclusions reached for 

this class will thus have some bearing on the problem of present interest. 

Subject to this restriction we have 

0 


(yr RU) ne 


(VV) =" (PRwv’). 
We shall study this relationship by considering simple examples, beginning 
with the case of constant radial flow: y R U =c. For this case the torque is | 


found to be =. 
rRu' v' =—crVv 


a 

& 1961 Energy Flows in a Vortex Tube ~ 353 
a 

on requiring that both velocity and mixing stress be zero on the axis. The 
work flux is 

4" 4*V Ruy =—cV?, 

inwards for c > 0, outwards for c < 0; that is, opposed in direction to the 


radial mass flow. The source term is 

a 

E 1 Te ey 
YY wt VRw u'v!) = y Or 


4 
ra 


Local heating occurs for inflow combined with swirl decreasing outwards and 
for outflow with swirl increasing outwards. The other combinations result in 

‘cooling. Practically, this means that radial inflow will result in the cooling of 
fluid in the core and the heating of that further out. 

: Having considered the effect of swirl on the work flux for a particular radial 
flow, we now take the other point of view. We take V = kr—*? and consider a 
variety of radial flows. Swirl distributions approximating to this form are 
commonly found in actual vortex flows; this was the case in the series of 
‘tests referred to here, and the agreement was particularly good for the con- 
figurations which exhibit the Ranque effect strongly. In any case, the quali- 
tative results obtained for this particular swirl distribution are more generally 
‘applicable; a specialized form is chosen for simplicity. We table the results 
for some simple examples: 


i Radial Flow Work Flux Source Term 
RU rV Ru'v' — = rv Rus) 
jcer =~ ta Y ae 
| (2) cr ord f 

(3) os ae x is 
B4)-c(a—r), rsa —ow (2-3) St se 
(5) cr(a—r), rsa : ce (f 7)? —oR(>- =) 


Case (1) models flow in the core; cases (2) and (3), that in the body of the 
flow; case (4), that near the tube wall. Case (5) covers the whole radius from 
axis to wall. Again we have chosen solutions such that the torque approaches 


zero in the core. 
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These examples follow the pattern that emerged on considering constant 
radial mass flow. Similar results are found for the more general swirl distribu- 
tion V = kr" where n > — 1. Stability considerations suggest that cases for 
which  < — 1 are not of physical interest. They were not found in experi- 
mental situations in the present study. 
We conclude that the radial work flux is usually in the direction opposed 
to the radial mass flux. But in the Ranque-Hilsch vortex tube the gas is” 
injected at the periphery and escapes through exits nearer the axis. Then 
radial flow must be predominantly inwards. We can expect that the core will 
commonly be cooled and the peripheral fluid heated through the working of 
the tangential stresses. 


9. Work Flux due to Axial Stresses 


The effect of the other important Reynolds’ stress is not simple. The 
approximate equation governing axial momentum has been shown [1] to be 


<2 Rw wv’) = eg Ley 


Requiring that the stress be bounded at the axis, we find the work flux to be . 


Y 


rW Roo! = —W fr dr. 
Oz 
0 


In order to demonstrate the properties that this function can be expected 
to possess typically, we consider a simple example: 


incorporating the property that W and 0P/dz are usually of opposite sign in 
the vortex flows studied. The work flux is 


rRW ww! = P(r —1) (r—3 


= 0 Tore Ly roe 
3 
<A (0) Siote a Cea 


This example illustrates the trends of more complex cases. Usually W and 


— / rv (OP/0z) dr have the same sign in an outer annulus, indicating an outwards 
U 


2 
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flux there. Most of the fluid in this outer annulus will be cooled; only that very 
near the wall will be heated. A complex flow pattern will be found in the core; 
it will not likely play an important part in overall energy separation. 


10. The Efficient Mechanisms of Energy Redistribution 


_ From the results of the four previous sections we conclude that the overall 
effect of each of the four important energy fluxes will be to transfer energy 
from the fluid in the core of the vortex to the fluid of the periphery. This 
agrees qualitatively with the behaviour of the Ranque-Hilsch vortex tube 
from which cooled fluid is withdrawn from the core at the inlet end of the 
tube. It was shown in Section 5 that the energy transfer mechanisms of 
turbulence are of the same order of magnitude as the observed energy redistri- 
butions. 

We have already seen that it is to the thermal flux Rw’ h’ and the work 
fluxes that we must look for important energy movements. The consideration 
of the details of the physical processes giving rise to these suggests the following 
conclusions regarding their relative importance. The efficient mechanism of 
energy redistribution in the core is thermal transfer by turbulent mixing 
across the large radial pressure gradients found there. The further separation 
‘in the outer fluid is primarily due to work fluxes associated with turbulent 
“mixing stresses. 

_ The analysis presented here gives a theoretical framework for the most 
plausible of the qualitative ideas on the Ranque effect that have been developed 
over the past thirty years. It is realized, of course, that this treatment contains 
such radical simplifying assumptions that it may accurately be described as 
‘crude. However, most of the assumptions have at least empirical support; 
many were, indeed, suggested by consideration of experimental results. It is 
felt that the conclusions reached should be applicable to the particular vortex 
tube configuration investigated by the present writer and probably for most 
of the configurations considered in previous studies of the Ranque effect. 


Symbols 
1,2 radial and axial co-ordinates ; 
u,v, W radial, tangential and axial velocity components (U, wu’, etc. are 


time-mean values and perturbations to them) ; 

0, ph, hy fluid density, pressure, enthalpy, and total enthalpy (R, 9’, etc. are 
time-mean values and perturbations to them) ; 

a, B, 0, €,7 quantities small with respect to unity. 
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Abstract ( . 


The energy separation within a vortex tube filled with turbulent compressible: 


fluid is investigated with an order-of-magnitude analysis of the energy equation. 
The physical processes corresponding to the important terms are: a heat flux: 
due to turbulent mixing of the compressible fluid through radial pressure and| 
temperature gradients, a flux of total energy produced by Archimedean} 
forces, and work fluxes associated with the two most important Reynolds’’ 
stresses. All these fluxes will commonly be outwards and will tend to cool | 
the vortex core. Experimental results are used to estimate the relative magni-- 
tudes of the contributions. The Archimedean éffect seems to be the least: 
important. 
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Zusammenfassung 


_ Die Energieverteilung in einem mit turbulenter kompressibler Fliissigkeit gefiill- 
ten Wirbelrohr wird durch Berechnung der Gliedergrésse der Energiegleichung 
untersucht. Die physikalischen Vorgange, die den wichtigen Gliedern entsprechen, 
sind folgende: ein Warmefluss, veranlasst durch turbulente Mischung der kom- 
pressiblen Fliissigkeit infolge der radialen Druck- und Temperaturgradienten; ein 
Fluss des Gesamtwarmeinhaltes, verursacht durch die Schwimmkraft; und die 
Arbeitsfliisse, die von den beiden wichtigsten Reynoldsschen scheinbaren Span- 
mungen abhangen. Alle diese Fliisse sind gewOdhnlich nach aussen gerichtet und 
dienen zum Kiihlen des Wirbelkerns. Die experimentellen Ergebnisse werden zu 
einer Abschatzung der Grdsse dieser Beitrage gebraucht. Der Schwimmkrafteffekt 
scheint am wenigsten wichtig. 


(Received: November 18, 1960.) 
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Note on the Large Deflection of a Circular Plate under a 
Concentrated Load 


By Sukumar Basutt, Calcutta, India’) 


Introduction 


For the large deflection of a plate we usually get non-linear equations which 
cannot be exactly solved. BERGER [1]”) has shown that if, in deriving the differential 
equations from strain energy, the strain energy due to second strain invariant in 
the middle plane of the plate is neglected, a simple fourth order differential equation 
coupled with a non-linear second order equation is obtained. He has solved these 
differential equations for the problem of circular plates under various boundary 
conditions subjected to normal uniform load throughout the plate. Though his 
approximation lacks sufficient technical interpretation for its justification, he has 
shown that his results for deflection, displacement, and stresses tally with the known 
results for all practical purposes. In this note following BERGER’s method, an 
attempt has been made to obtain the deflection for a circular plate under a con- 
centvated load at the centre. 

1. In the case of circular symmetry if / is the thickness of the plate, w the 
displacement perpendicular to its middle plane, « the radial displacement in the 


1) Brahmananda Keshab Chandra College. 
2) Numbers in brackets refer to References, page 362. 


Se 


358 Kurze Mitteilungen — Brief Reports — Communications bréves ZA 
middle plane under concentrated load at the centre, then, following Bence 
method, the differential equations for w and u will be 


a 1 d\(dw 1 dw 4, )_ | 
(ars Y ) a , oo w) =0, (1.1) 


Aen ye ae 


except at the load and | 
du U 1 (2 DOT apt a 2y 
dy Y 2\ar ) ¥ i 


ie 


where « is a constant of integration. 

Again considering the radial stress and shearing stress on a concentric circular 
ring of radius y and the concentrated load P at the centre we have [2], as v tends 
to zero 


ad (dw 1 dw Eh du 1 /dw\2 u| dw P 
Ds (Fa Y 7) = Ta le+ale) te |a+ae 1.9 
where D = Flexural rigidity of the plate 
= E h'/12 (1 — 0), 
o = Polsson’s ratio, 
E = Youne’s modulus. 
Substituting (1.2) into equation (1.3), we have the relation 
ad (@w 1 dw Eh dw J 
re: a 
D1, (Fa Y ar a? w) fg era 
as y tends to zero. 
At the centre of the plate both u and dw/dr will be zero; hence we have 
‘ ad (dw 1 dw P 
= 22 od x 
pet oad dy ( dr® a ee w) 2% (C7 


2. A circular plate under concentrated load at the centre with clamped edge.. 


Let the radius of the plate be a. Then fora clamped edge the boundary conditions; 
will be 


dw 
Basan (0, u=0O when rv=a. (2.1) 


Let us take w in the form 


w=AT,(ar) +B (Kolo ”) + log) 4e (2. 2) 


where I(x 7) and K,(« 7) are the modified Bessel functions of the first and second 
kind of order zero. Substituting this w in (1.4) we have 


we 
B Se On D ea . (2: 3) 
Using boundary conditions for w and the relations [3] 
Iy(ar) = Iy(ar), Ké(ar) = —K,(« 7) 
and 
1 


Iy(ar) Ko(a7) + Ig(a 7) Ky(a 7) = = 


: 
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we get 

= K,(aa) —1 Iy(a a) —1 

Z Ae B Bec HES 

; Serene Co PAS lhaaddoa, ee) 
4 faking these values of 4, B and C, w is given by 

; P 

4 aay a Fis “aa[K,(« a) Ip(a7v) + Ky(a 7) I,(« a)] 

+ aa I,(a a) log” — Io(ar) + Io(n a) — 1}. (2. 5) 
J To find the displacement u we have from (1.2) and (2. 2) 

. 2 je . 

; ur= = yr 7 | [4287 Ben + Bit Kiar) 

: —~2A Bo®rI,(ar) Kiar) +24 Ba I,(x7) 

E , — 2 Ba K,(ar i+ ]a+F 

‘where F is a constant of integration. 

Integrating these expressions involving Bessel functions [4] 

: 2 2 2 

ny = cae v2 — So tz a2 v2 [I2(a7) — R(ar)] + ar Iy(a7) Ip(x n| 

4 24 2 2 

, B? (1 

: ee meas es a? v? [K2(a 7) — Ki(a7)] — ar Ky(a 7) Ko(x n| 

1 

‘ +AB fe a2 v2 [T,(a7) Ky(a7) + Ip(a 7) K(x 7)] | (2. 6) 


— ar (Ly(ar) Kolar) — Tolar) Kile m1} 


2 
—ABI,(a7r) — B? K,(«7) — rote logy + F 


“where B and A are given by (2.3) and (2.4). 
Using the boundary condition for u 


; B2 (1 1 — I,(« a) anh (a4) a? fh? a 
, => {> —s = pea? 
ee ey 3 SHOE a Tiaa) 2 \ Tyla 2) 24 end} 
To determine the constant «, we shall use the condition that u = 0 as 7 > 0, which 
gives : ih 2% 
| B ib o aa K,(a a) — 
re Jee ae F=0 
2 2 + ¥ + log, aa I,(« a) \ 4 
or 
a! 
ee 6 

Pa )’ x 3 ea} ~ (2.8) 
ag Ee es I,(aa)+aakK,(xa)—2 1 (Zale2) — =)" 
; Peeeh 2 a @ Ty (« a) 2\ Iya a) 


where y is EULER’s constant. 
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The deflection will be maximum at the centre of the plate and it is given by 
P aa Ne a | 


i 
| 


= _ Pap ) 
290 Dot" {v + log Z aa I(x a) ( 9)) 


As « tends to zero the equations (1.1) and (1. 4) reduce to the corresponding equations: 
for the problem of small deflections of circular plates under concentrated load at: 
the centre. 

Taking the limit of (2.5) as « tends to zero we have [5] 


: | 
= es, ( logy 13 (a*— 1), (2.10)) 


CS 


which is the expression for small deflection of a circular plate under concentrated! 
load at the centre, the boundary being clamped. 

For a thin plate flexural rigidity of the plate will tend to zero and « will be: 
infinitely large. Expanding the Bessel functions asymptotically [6] in (2.8) and. 
(2.9) and taking the limit as « tends to oo, we have on the assumption that h? «2! 
exists, the central deflection for a thin membrane represented by the term 


Li Wo sn ( 1— oe 1/3 
aro (Pa@jEhi® —\2nx(y— 1/2 ee 
Putting o = 0-3 we get 

Ht" (Pale His 


In the case of a membrane the approximate formula deduced by TIMOSHENKO, 
Wotnowsky-KRri&GER [7] leads for the deflection at the centre to the result 


O-713% (2.11). 


Wo 
(P @/E hj 


= 0-789. (2.12) | 


For values of « other than the extreme ones discussed above, a graph (Figure 1) 
has been drawn showing w,/hk against P a2/x D h. 


3. A circular plate under concentrated load at the centre with supported edge. . 
At the supported edge deflection, the bending moment and displacement normal 
to the boundary will be zero. In this case the boundary conditions will be 


dw d 
as Ae te, “w#=0, when v=@a@. (3.1))) 


Taking the same expression for w as before and using these boundary conditions, 
the arbitrary constants are given by 


P 
See igaieie (Sg 
H K,(« a) + aaK,(«%a)—1 
Aes es 
| aa ode Mae ae Ca 
thr 
aa Tua) Maite Qior (3.4) 
where 
a? a2 


Ee eee (3.5) 
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ee eh a 


Substituting these values of A, B and C in the expression for w and after 
‘simplification we have 


: P 
e 2nD|H I,aa)—aal,aa)] 2 aUGleahistee) goin e) tale a] 


+ aa I,(a a) log” Seelslea?) oe Lalani | (3.6) 


+H [Kole a) I,(a 7) — Ky(a7) I(x a) — Iq(« a) log~ |! . 


45 
40 LINEAR THEORY 
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35 
KNOWN APPROXIMATE SOLUTION 
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ee 
She 
2-0 
V5 PROPOSED SOLUTION 
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SUPPORTED FOR © =0:3 
7-0 
0-5 
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Figure 1 
Values of a showing w,/h against P a®/a Dh. 


The displacement w in the plane of the plate will be given by the expression (2.6) 
where as B, A and C are given by (3.2), (3.3) and (3.4) respectively. 
Proceeding exactly as before 


(H_— 1) (1o(« 4) — 1) 
H I,(a a) — «a I,(a a) 


a [H Ia a) —aaT,(xa)) 


i 
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and the equation to determine « will be 


= (xa) 
eaali= T,(« 4) Kies) 2—H[I,(aa)—K(wa)—1] || 
Diet) (aa) +aa “a)—2— o(a a) — = at 
i ya Tey ere Higgs 4 
H2 oa? { 

a [z+ fy see | [Io(« a) — 1]? 4 

; [H I)(a a) — aa I(x a)]? : . 


The central deflection will be given by 


P a a Tole d) 406 Bilao) +2 Bola) — 2) 7 | 

wo 2% D ot {7 + log a HI Iy(« a) — «a I,(« a) ‘ (3.9)} 
Again as « tends to zero equation (3.6) reduces to 

ee SFEAO 7 ra - “] 

= 2 rv? log — 3.10) | 

16 aD rast Ser ee cae | 


a well known result for small deflection. As « tends to infinity the same result: 
(2.11) can also be obtained from (3.8) and (3.9) to give the central deflection for’ 
a membrane with supported edge and loaded at the centre. In the same figure: 
referred to before w,/h has been plotted against P a*/x Dh with o = 0-3. . 

In conclusion, the author wishes to thank Prof. B. SEN, D.Sc., F. N: agi 
Jadavpur University for his guidance in the preparation of this note. 
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Zusammenfassung 


Es werden die grossen Deformationen einer Kreisplatte unter zentrischer Ein- 
zellast nach der Methode von BERGER ermittelt. 


(Received: 11 November, 1960.) 
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Non-Steady Periodic Boundary Layer 
iy By Miss DurGa Roy, Calcutta, India 


Introduction 


_ For many years people are trying to study the periodic phenomenon of the wake. 
4 It is actually the periodicity of the wake for which it is needed to study the 
non-steady character of the motion. Calculations on non-steady boundary layer 
equations have been performed by many authors. Biastus [1]*) in 1908 first showed 
the approximation method of solving non-steady boundary layer equations. In 
the same year Bottze [2] performed the calculation on steady boundary layer 
equations of an axi-symmetric body of revolution. It is however also very in- 
structive to study the development of boundary layer periodic in time. In the year 
1932, SCHLICHTING [3] found out a solution periodic in time by the approximation 
‘method of the two dimensional non-steady boundary layer equations. In this paper 
the same method is applied to find out a solution periodic in time of non-steady 
boundary-layer equations of an axi-symmetric body of revolution. The only dif- 
ference is that in this case a different series expansion for the stream function is 
taken. 
. Boundary Layer Equations and its Solution 
It is shown in the case of two dimensional periodic boundary layer that the 
‘successive approximation method may be used in cases when the amplitude of 
‘oscillation is small compared with the dimensions of the body. This condition is 
‘satisfied in the two dimensional case. It can be easily verified that this is satisfied 
in the axis-symmetric case too. 

The boundary layer equations of an axially-symmetric body of revolution are 
given by, 


Ou i 
Ot 


where (x) is the radius of the cross-section of the body of revolution. 
The equations (1) are to be solved with the boundary conditions, 


Ou Ou 1 Op i Pu 0 


0 
= - —(vv) = 0. 1 
Ox ne Oy Q Ox 4 Oya Ox ee (1) 


oy 


y=0, “#=0=0, Y= O, Maa OA t) (2) 


where U(v, t) denotes the potential flow about the body of revolution. 
The pressure impressed on the body follows from the non-steady BERNOULLI 


equation, 


(3) 
Applying the general approximation method (1) infinding out a solution of the 

non-steady boundary layer equations we can write 
u(x, y, t) = U(x, y, t) + uy(%, ¥, 2) - (4) 


1) Number in brackets refer to References, page 366. 
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Thus from (4) and (1) we obtain the first and second approximation equations 
which are given by | 
OUy 079 0U Ou, 07, 0U OUy OUy | 

ae = Ey ff peels amt yay Satis Ney BN BN ee 5) 

be a a8 RT eee ic de emily ris ek 

j 

with boundary conditions ; 
y= O, t= t= 0, yeoo, w= U4, u4=0. (6) | 


Let us take the potential flow to be periodic in time, i.e. 
Cat) Ua) ee 


We define 7 = y //n/v a dimensionless quantity. 
Let 


on )= //2 te x) 7 a(n) emt (7) 


: U, ne © 
int v 0 nt 
é =I /2 * Coe? a (8) \ 


x 


hence, 


Uy = U, Ch e?2", 


From (8), (5), and (6) we get _ : 
we too— Co = 2 (9) 
with boundary condition 
=O at w=0, QO=1 at n= co. 


The solution of (9) is given by 


fit apel ae (10) 
Taking the real part we get . 


uy =U [cost — exp(— 1) cos(n t alk (11) 


For the second approximation let the stream function be of the form 


wiles ot) = |/% [oy 220 ce atny etme + Saale 
+ ae Cielo ett + Calms]. (42) 


n ax 


From (7), (8), and (12) we write the real expression for Up, Vo, and %y, i.e. 


ea a Pet n ay OU, U, or | 
My = (Coet™et Coe int), vy = hee ert - we = = ae (13a) 
= iy 0U 1 re 2int ¥ —Qint / yr 
eee > os a (Sia é + lige T Cup fF Gint 
ON 2int ae e-2int ' 
Wp tap ALi are of ab tiet cial ee 


Using (13), (5), and (6) we get 


rary wy 1 , " 
So Shy Sid (1 — Oo? + a Co) » (14) 
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= 24 Ci¢ te=— ae , (15) 
: it: puned tar Sabea Une ok 

: — ay= er Co So + 4 (Gg Co Gy Go)s (16) 
4 1 

4 — Sta = — > (So So + So So) - (17) 


The normal and tangential components of the periodic contribution must vanish 
at the wall, whereas at a large distance from it only the tangential component 
vanishes. 

__ Equations (14) and (16) are exactly the same as in the two dimensional case so 
the solutions are known. Equations (15) and (17) are solved actually. Putting 


a = nlV2 


Bie = 1 expt 1+.) 11 - FS expl- 0 +) 11 


: — J exp[— 29 (L4+9)+ Gi expl-y2x' (1+ a). (18) 
tia= F + en’ — 1) — 6 (y+ Qexpl- (0 +4) 71 
bs 5 [(n’ — 1) +4 (n’ + 2)] exp[— (1—4) 71], (19) 


if 
Ge=Cia=9 at n=0, tic= 9, fia. mie at n=O. 


It is found that only the boundary conditions at the wall can be satisfied, and 
that at a large distance the tangential component is finite but not zero. 
_ Thus from the second approximation we see that at a large distance from the 
body i. e. outside the boundary layer the velocity does not vanish. 
Its magnitude is given by 
3 OU, Cena, 


iigta, OS, 9) = 4n Uo Ox 2nv dx° 


(20) 


So we see the remarkable result that a potential flow which is periodic with 
respect to time induces a steady secondary motion at a large distance from the 
wall of the body of revolution. This steady secondary wz is composed of two terms 
and this u, is the same as in the two dimensional case if we neglect the second 


‘term of it. 
Its magnitude given by equation (20) is independent of the co-efficient of 
‘viscosity but depends on the shape of the body, i.e. the radius of the cross-section 


of the body of revolution. 
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Zusammenfassung | 


Es wird die nichtstationare Grenzschichtbewegung an einem Rotationskérper 
untersucht und die Stérungsgeschwindigkeit bestimmt. Ahnlich wie im ebenen Fall | 
verschwindet diese in grdésserer Entfernung nicht und ist stationar. 
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Die Fernrohre und Entfernungsmesser. Von A. K6nic und H. K6HLER 
(Springer-Verlag, Berlin 1959, 3. Auflage). 475 S., 471 Fig.; DM 88.-. 

Es handelt sich um eine durch H. Kéhler vollstandig neu bearbeitete Auflage 
des bekannten Buches von A. Konic. Der erste und umfangreichste Teil behandelt — 
zunachst die theoretischen Grundlagen und dann die wichtigsten heute gebrauch- 
lichsten optischen Systeme, wie Hand-, Beobachtungs-, Zielfernrohre, geodatische, 
feinmesstechnische, astronomische und photographische Fernrohre und enthalt 
ausfihrliche Tabellen tiber Linsenobjektive, aplanatische Spiegelsysteme und Oku- 
lare, ebenso wurden mit einer gewissen Vollstandigkeit die verschiedenen Priif- 
methoden aufgenommen. Nach einem kurzen zweiten Teil iiber Mikrometer, der 
praktisch unverdndert aus den friiheren Auflagen iibernommen wurde, folgen noch | 
die verschiedenen Typen von optischen Entfernungsmessern, wie sie fiir militarische, 
geodatische und andere zivile Zwecke Verwendung finden oder fanden, Anhangs- 
weise finden sich noch Zusammenstellungen der wichtigsten Beziehungen der 
GauBschen Dioptrik und der lichttechnischen Einheiten, ebenso ein Abriss der 
Geschichte des Fernrohres und Literaturangaben. Das Buch, das weder zu sehr auf 
technische Einzelheiten noch auf theoretische Entwicklungen eingeht, stellt eine 
ausgezeichnete und empfehlenswerte Einfiihrung in die praktische Optik dar. 

E. RotH-DESMEULES 


Sammelband zur statistischen Theorie der Turbulenz. Mit Beitragen von 
A. M. Kotmocorovy, A. M. Osucnow, A. M. JAGLom und A. S. Montn. In deutsche 
Sprache iibertragen von Dr. H. GorRING (Akademischer Verlag, Berlin 1958). 
228 5., 27 Abb.; DM 25.50. 

Die statistische Darstellung der Turbulenz der Strémungsfelder von Fliissig- 
keiten und Gasen ist in den letzten 15 Jahren durch die grundlegenden Beitrage 


von G. T. TAYLor, von Karman, HEISENBERG, KoLMOGOROV und anderen weit- : 


gehend entwickelt und geférdert worden, Ein massgeblicher Anteil kommt den 
russischen Arbeiten zu. Da iiber diese Ver6ffentlichungen in den bis heute erschie- 
nenen zusammenfassenden Abhandlungen und Lehrbiichern oft nur sehr sparlich 


i i 
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berichtet wird, ist es lebhaft zu begriissen, dass der Akademische Verlag, Berlin, 
einen Sammelband tiber die bedeutendsten russischen Arbeiten herausgegeben hat. 
‘Die Ubersetzung ist in klarer Formulierung durch Dr. Herbert Goering ausgefiihrt 
worden. 
| Unter den wiedergegebenen Publikationen befinden sich die grundlegenden 
_Arbeiten von A. M. OsucHow iiber die statistische Beschreibung und Struktur 
turbulenter Felder. Ferner die Notizen von Kotmocorov iiber die Entartung 
isotroper Turbulenz in einer inkompressiblen Fliissigkeit sowie die Zerstaubung 
-von Tropfen in einer turbulenten Str6mung. Schliesslich gebtihrt der Arbeit von 
_A. S. Monin und A. M. OBpucuow iiber die fundamentalen Gesetzmassigkeiten der 
turbulenten Vermischung in der bodennahen Atmosphare besondere Erwahnung. 
R. SANGER 


Wirbelstré6me und Schirmung in der Nachrichtentechnik. Von H. Ka- 
DEN (Technische Physik in Einzeldarstellungen, Bd. 10) (Springer Verlag, Berlin 
1959). 354 S., 195 Abb.; DM 66.-. 

Das vorliegende Buch ist eine vollstandig umgearbeitete und erganzte Neuauf- 
lage der 1950 erschienenen Monographie «Die elektromagnetische Schirmung in der 
-Fernmelde- und Hochfrequenztechnik». Alle Untersuchungen sind fir stationare 
-Wechselfelder durchgefiihrt. Fussend auf den Maxwellschen Gleichungen in Vektor- 
-schreibweise und unter Benutzung angepasster Koordinatensysteme, werden in 
zwolf Kapiteln etwa folgende Probleme durchgerechnet und Naherungsverfahren 
dazu abgeleitet: Wirbelstréme in Blechen, Drahten und Kugeln; Schirmwirkung 
-metallischer Hiillen gegen dussere magnetische Wechselfelder in Form von ZWEl 
parallelen Platten, Hohlzylinder und diinnwandiger Kugel aus magnetischem und 
“unmagnetischem Material; Riick- und Schirmwirkung metallischer Hiillen gegen 
linear polarisierte elektromagnetische Wellen (Radar); metallische Hiillen mit 
innerer Felderregung, zum Beispiel koaxiale Leitung, Doppelleitung und Stern- 
-vyierer mit koaxialem Metallmantel, Spulen in Metallhiillen; mehrschichtige Schirme 
aus verschiedenen Metallen; Durchgriff von elektrischen und magnetischen Feldern 
durch Spalte, Schlitze und Locher, was zu Kopplungen zweier Leitungen benutzt 
wird (Richtkoppler der Mikrowellentechnik) ; Wirkung begrenzter offener Schirme; 
Schirmwirkung von Drahtgittern; Schirmung gegen Stdérstréme bei Durchfih- 
‘rungskondensatoren, Doppelkontaktbuchsen und -federn; Kopplungswiderstand 
von Leitungsschirmen (wichtig fiir das Nebensprechen). Das Schlusskapitel orien- 

‘tiert iiber wichtige Eigenschaften der Zylinder- und Kugelfunktionen. 
So reichhaltig auch der Inhalt erscheint, es hat alles einen starken innern Zu- 
sammenhang, was das Studium einer Einzelfrage sehr erleichtert. Das Buch dient 
‘aber nicht nur dem Theoretiker, sondern durch seine zahlreichen numerischen 
Diagramme und Tabellen ebenso dem praktizierenden Ingenieur, der mit Schir- 
“mungsaufgaben zu tun hat. Dem Buch ist eine recht weite Verbreitung zu wiin- 
-schen, denn in seiner Vollstandigkeit und Abgeschlossenheit steht es wohl einzig da. 
H. WEBER 


Report on the Algorithmic Language ALGOL 60. Von PETER NAuR 
(Acta Polytechnica, Copenhagen 1960) 40 S.; Sw.-Kr. 7.-. 

Die algorithmische Sprache ALGOL wurde geschaffen, um Rechenprozesse 
der numerischen Mathematik in einer Form zu beschreiben, welche einerseits 
_méglichst lesbar sein sollte, andererseits aber die automatische Ubersetzung in ein 
-Maschinenprogramm durch den Rechenautomaten selber erlaubt. Der vorliegende 
‘Bericht enthalt, ausser einigen einfiihrenden Bemerkungen, die definierenden Re- 
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geln der Sprache in konzentrierter Form uid ist keineswegs als Einfithrung fiir 
den Beniitzer gedacht. Es handelt sich im iibrigen um einen Abdruck aus Nume- 
rische Mathematik 2, 106 (1960). . P. LAUCHLI 


VIII. Internationaler Astronautischer _Kongress - Barcelona 1957. 
Berichte — Proceedings — Comptes rendus. Herausgegeben von P. J. BERGERON, 
W. von Braun, K. A. Exricke, F. Hecut, J. MARIAL, E. SANGER, K. SCHUTTE, 
L. I. SEpov, L. R. SHEPHERD, S. F. SincErR. Schriftleitung F. HeEcur (Springer- 
Verlag, Wien 1958). VII + 607 S., 233 Fig., 91 Tab.; SFr. 128.-. 

Der VIII. Kongress der Internationalen Astronautischen Féderation (Barcelona, © 
6.-12. Oktober 1957) fand statt unmittelbar nach dem Start des ersten Erdsatelliten | 
«Sputnik I» (4. Oktober 1957), dem erfolgreiche Abschiisse weiterer Satelliten 
folgten. Es war sowohl technisch wie geschichtlich eine sehr interessante Zeit- 
periode von angestrengten Bemiihungen um den effektiven Beginn der Raumfahrt 
im Rahmen des Internationalen Geophysikalischen Jahres 1957/58. 

Einen interessanten Einblick in den Stand der damaligen Arbeiten vermittelt 
der vorliegende Bericht. Er enthalt 46 Beitrage in deutscher, englischer, franz6- | 
sischer, italienischer und russischer Sprache, jeweils begleitet von anderssprachigen 
Zusammenfassungen. Die behandelten Themata lassen sich etwa folgendermassen | 
einteilen : . 

I. Allgemeines tiber Astronautik — Referate tiber wissenschaftliche und tech- 
nische Ziele der Raumfahrt, wie u.a. Erforschung von Sonne, Planeten 
sowie von intelligentem Leben auf Planeten; 

IT. Erforschung der hohen Atmosphare — Resultate der Atmospharenmessungen - 
durch Ballone und Ballontechnik; 

III. Raumfahrt- und Astro-Physik — Schutzmassnahmen gegen Andruck, Dyna- 
mik der Raumk6rperbewegungen und die zugehérigen Berechnungsmethoden 
sowie Programmierungen, relativistische Raketenmechanik, Meteoriten, | 
Fragen hoher Temperaturen sowie Untersuchungen von Raumen mit Da- 
seinsbedingungen fiir belebte Organismen (Okospharen) ; . 

IV. Raumfahrttechnik — Probleme der Riickkehr von Satelliten, Steuerung, © 
Stufenprinzip, Lastenoptimalisierung, Projekt einer dreistufigen bemannten 
Kleinstrakete «Meteor- Junior»; 

V. Raketenantriebe — Stabilitat der Fliissigkeitsraketen-Motoren, Variation des 
Raketenschubes, Stré6mungen von Zweiphasengemischen, Nuklear- und 
Ionenantriebe ; 

VI. Kiinstliche Satelliten — Studien iiber Gewichtsaufwand fiir Satelliten, An- 
wendungsméglichkeiten von Satelliten zu Kommunikations- und Messzwek- | 
ken sowie Feststellung und Beobachtung von Satelliten selbst; 

VII. Raumfahrtbiologie und -medizin — Herstellung des gewichtslosen Zustandes 
in Fliigen, biospharischer Index, Einfliisse einer totalen Sonnenfinsternis 
aufs Blut; 

VIII. Raumfahrtrecht — Entstehung von neuen Rechtsfragen durch die Raumfahrt 
bedingt. | 

Die Vielfalt und der Umfang des Dargebotenen verunméglichen es, hier auf die 
Einzelheiten einzugehen. Der interessierte Leser findet im Barcelona-Bericht reich- 
liches Studienmaterial, das einen zweckmassigen Anschluss an die inzwischen 
weiter vertieften bzw. neugewonnenen Erkenntnisse und Erfahrungen in Astro- 
nautik samt verwandten Gebieten gestattet. Eine sehr gediegene und iibersichtliche — 
Gestaltung des Buches erleichtert das Studium wie das Nachschlagen der einzelnen 
Beitrage. Z. PLASKOWSKI 

: 
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Eine Einfiihrung in die Eigenschaften . Band 19 


_und Probleme der Kernreaktorwerkstoffe - Sammlung «Lehr- und Handbicher 
‘ der Ingenieurwissenschaften» 


407 Seiten mit 141 Figuren 
(1961) Ganzleinen Fr./DM 57.50 


‘Inhalt: » 


. Einfiihrung in die Reaktortechnologie 

, Physikalische Grundlagen; Aufbau und Einteilung wichtiger Kernreaktortypen 
_ 2, Fiir den Reaktorbau wichtige Eigenschaften der Materialien aoe 
_ Atomeigenschaften “hy 
Struktur- und Aggregatseigenschaften von kristallinen Festkérpern 

Eigenschaften der Fllssigkeiten und Gase 

£ Oberflachenreaktionen von Festkérpern mit Gasen und Flussigkeiten ‘ : d 
_ 3. Die Wirkung der Reaktorstrahlungen auf die Materialien ' r 
4 Natur, Energie und Intensitat der in Kernreaktoren auftretenden Strahlungsarten 

- UWberblick Uber die Wechselwirkungen zwischen Strahlen und Festkorpern : 
E Die Stérung des Kristallgitters durch Bestrahlung im Kernreaktor und damit verbundene Eigenschafts- ae 
__veranderungen 4 tes 3. 
if Erwarmung von Festkérpern durch Strahlen in Kernreaktoren ‘ ne ; 
4 “a 
4 

; 


Chemische Wirkungen der Strahlungen in Reaktormaterialien ‘ 
Die induzierte Radioaktivitat in Reaktormaterialien 7 : 
Die Spaltprodukte 


_ 4, Obersicht tiber die Reaktorwerkstoffe Lys en! 
Spalt- und Brutstoffe; Hillstoffe flr feste Spaltstoffelemente; Moderator- und Reflektormaterialien; te 


: KuhImittel; Absorptionsmaterialien; Konstruktionsmaterialien; Abschirmmaterialien 


_ Im vorliegenden Buche wird nach einer kurzen Einleitung in die Technologie der Kernreaktoren eine Ein- 

- fihrung in die moderne Theorie der Werkstoffeigenschaften gegeben. Dabei werden im Kapitel 2 anhand 

von in der Kerntechnik wichtigen Beispielen die Eigenschaften der Atomkerne, der Atomhiille, der Kri- 

stalle und vielkristallinen Stoffe besprochen und das flr die Technik interessante Verhalten der Metalle, 
Legierungen, Fliissigkeiten und Gase beleuchtet. Im dritten Kapitel wird die Beeinftussung der Festkorper nod 
durch energiereiche Strahlen behandelt, wobei wiederum speziell jene Beispiele hervorgehoben werden, ; b 
~ welche fiir den Reaktorbau von Interesse sind. Das Schlusskapitel gibt einen Uberblick Uber die heute 1 
_ gebrauchlichen Werkstoffe und KiihImittel von Reaktorkernen. Es wurde absichtlich darauf verzichtet, es 
die Herstellung von Werkstoffen und die Fabrikation wichtiger Konstruktionsteile, zum Beispiel der Spalt- { 
stoffelemente, zu behandeln. Das vorliegende Buch soll vielmehr eine griindliche Einarbeit in die Werk- 
stoffprobleme der Kerntechnik erméglichen und-die Stoffeigenschaften sowie deren Ursache und Ver- 

anderung unter Bestrahlung.nach dem Stande der heutigen Kenntnisse erhellen und beschreiben. Es 

wendet sich vor allem an Ingenieure, daneben wird jedoch mancher Abschnitt auch Chemiker, Physiker 

und Metallurgen interessieren, welche auf dem Gebiete der Kerntechnik tatig sind oder sich mit den in 


ihr auftretenden Problemen vertraut machen méchten. 
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' Neuerscheinung — 


Statistische Methoden 


Aaa = 
fiir Naturwissenschafter, Mediziner und Ingenieure 
‘ von 


- ARTHUR LINDER 
Professor an der Universitat Genf und an der Eidgenéssischen Technischen Hochschule Ziirich — 


Dritte, umgearbeitete und stark erweiterte Ausgabe » 
' (1960) 484 Seiten mit 58 Figuren und 73 Beispielen. In Ganzleinen Fr. 54.— (DM 54.-) 


Mathematische Reihe, Band 3, Sammlung «Lehrbiicher und Monographien aus dem Gebiete der exakten Be 
as s ’ Wissenschaften» ‘ : 


Angesichts der unvermindert anhaltenden Nachfrage hat sich der Verfasser entschlossen, seine bekannte 
. Einfiihrung in die neueren Methoden der mathematischen Statistik stark zu erweitern. Gleichzeitig. wurde 
' der den Anwendungen gewidmete Teil des Buches vollig umgestaltet. i 
Neu aufgenommen wurden unter anderem: Verschiedene Anwendungen von Chiquadrat, die Streuungs- 
zerlegung bei ungleichen Klassenhaufigkeiten, die Bestimmung von Streuungskomponenten, die nichtlineare 
Regression, die «analysis of covariance», das Schatzen von Parametern und einige Transformationen von 
Prozentzahlen (Arc sin, Probit, Logit, Loglog). . 


1959 erschien vom gleichen Autor: 


Planen und Auswerten von Versuchen > 


Eine Einfiihrung fiir Naturwissenschafter, Mediziner und Ingenieure 
(1959) 2. Auflage, 182 Seiten mit 9 Fig., Ganzleinen Fr./DM 21.- 


Reihe der experimentellen Biologie Band 13, Sammlung «Lehrbiicher und Monographien aus dem Gebiete 
: der exakten Wissenschaften» : 


Diese Einfiihrung richtet sich an Naturwissenschafter, Mediziner und Ingenieure; sie setzt keine Kenntnis 

_ in mathematischer Statistik yvoraus. Der Leser wird angeleitet, Versuche richtig zu planen und einwandfrei _ 
auszuwerten. Die griindlich durchgearbeiteten Anwendungsbeispiele aus der biologischen, medizinischen, 

_ industriellen und landwirtschaftlichen Forschung bilden einen Wichtigen Bestandteil des Buches. 


Kiirslich erschien vom gleichen Autor: 


Handliche Sammlung 
mathematisch-statistischer Tafeln 


(1961) 40 Seiten, Fr./DM 4.50 


Die yorliegende Sammlung statistischer Tafeln soll die umfangreicheren, 
‘sie stellt vielmehr einen handlichen Auszug der am meisten verwendeten T: 
stets griffbereit auf dem Arbeitstisch zu liegen. Ausser den Tafeln von ¢ 
Tafeln fir die Transformationen von Prozentzahlen mittels Arcus sinus 
_ Tafeln fiir zufallig angeordnete Zahlen, 
zu beniitzen sind, 


bisher erschienen nicht ersetzen: 
afeln dar, die dazu bestimmt sind, 
» x* und F enthalt diese Sammlung 
, Probits, Logits und Loglog, sowie 
die besonders beim Planen von Versuchen und Stichprobenerhebungen 


z. Zu beziehen durch Ihre Buchhandlung — Obtainable from your bookseller 
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